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CHAPTER  1 


INTRODUCTION 


1.1  Background 

1.1.1  History  of  Auroral  Science 

An  excellent  description  of  the  history  of  auroral  science  was  given  in 
a  lecture  by  the  great  Sydney  Chapman  at  the  NATO  Advanced  Study 
Institute,  held  at  As,  Norway,  July  29  through  August  9,  1968.  This 
summary  was  derived  from  the  publication  of  Chapman’s  and  other 
lectures  put  together  by  McCormac  and  Omholt,  in  Atmospheric 
Emissions. 

Chapman  supposed  the  beginnings  of  auroral  history  dated  back  to 
1600,  when  Gilbert  discovered  the  earth  is  a  great  magnet.  In  1621  a 
magnificent  auroral  display  was  witnessed  by  Gassendi  from  the  south  of 
France;  he  was  perhaps  the  first  scientist  to  view  it  objectively.  He  named  it 
the  aurora  borealis,  or  northern  dawn.  Nearly  a  century  later,  in  1716, 
Halley  was  the  first  to  suggest  a  connection  between  the  aurora  and  the 
geomagnetic  field,  when  he  observed  the  great  aurora  in  London. 
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Also  in  London,  in  1722,  Graham  discovered  the  transient  variations 
of  the  earth's  magnetic  field,  which  are  regular  on  some  days  and 
disturbed  on  others.  Celsius,  with  an  instrument  from  Graham,  began 
similar  studies  in  Uppsala,  Sweden.  By  correspondence  they  found  that 
days  disturbed  at  Uppsala  were  disturbed  also  in  London.  In  1741  Celsius 
and  his  assistant  Hiorter  noted  that  the  occurrence  of  aurora  is  linked  with 
magnetic  disturbances.  Thus,  when  Gassendi  in  France  and  Halley  in 
London  saw  an  aurora,  there  must  have  been,  unknown  to  them,  what 
Humboldt  in  1806  called  a  magnetic  storm.  A  reasonable  estimate  of  the 
height  of  an  aurora  seen  in  England  in  1784  was  made  by  Cavendish. 

In  1833  the  geographer  Mucke  recognized  the  existence  of  the  auroral 
zone  of  maximum  frequency  of  auroral  visibility,  from  his  studies  of  the 
auroral  reports  of  arctic  explorers,  and  from  the  long  records  made  by 
atomic  chemist,  Dalton.  By  1844,  Schwabe  had  recognized  the  occurrence  of 
the  sunspot  cycle,  and  within  a  decade  it  was  found  that  the  magnetic 
variations,  both  regular  and  enhanced,  are  more  frequent  around  sunspot 
maximum.  Owing  to  the  poor  quality  of  auroral  statistics  it  took  decades 
before  it  was  established  that  the  aurora  also  follows  the  sunspot  cycle. 

In  1846  the  Scottish  geophysicist  Broun  thought  horizontal  magnetic 
measurements  had  a  lunar  periodicity,  but  in  1858  he  re-interpretted  the 

o 

variation  as  being  connected  with  the  solar  rotation.  In  1868,  Angstrom 
made  the  first  observations  of  the  auroral  spectrum,  the  beginning  of  a 
great  chapter  in  auroral  science,  but  it’s  study  made  little  progress  in  the 
nineteenth  century  because  of  what  atomic  physicists  were  yet  to  discover. 
In  1882,  the  Scottish  physicist,  Stewart,  inferred  the  existence  of  what  we 
now  call  the  ionosphere,  concluding  that  the  dependence  of  the  daily 
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magnetic  variation  on  the  sunspot  cycle  could  only  be  explained  by  electric 
currents  in  a  conducting  layer  of  the  upper  atmosphere. 

The  advances  toward  an  auroral  theory  awaited  the  discovery  of  the 
electron  by  atomic  physicists,  around  1880,  and  also  on  experiments  and 
theory  that  showed  how  the  electron's  motion  is  affected  by  magnetic,  as 
well  as  electric  fields.  This  led  Birkeland  in  1896  to  infer  that  auroras  occur 
in  polar  regions  because  the  geomagnetic  field  must  deflect  migratory 
electrons  toward  those  regions.  This  and  other  theories  were  important 
beginnings  in  auroral  science,  and  were  greatly  modified  through  the 
years. 

One  of  the  most  notable  areas  of  auroral  study  which  has  made  great 
progress  is  in  spectroscopy  of  the  aurora.  Measurement  and  identification 
of  numerous  atomic  lines  and  molecular  bands  in  the  auroral  spectrum 
was  a  great  task  undertaken  by  many  scientists,  and  is  still  ongoing. 
Though  many  unknowns  still  exist,  a  great  deal  is  known  about  the  aurora, 
thanks  to  studies  and  advances  in  instrumentation  during  the  20th 
century. 

1.1.2  Fundamentals  of  Aurora 

Auroral  activity,  as  implied  above,  has  a  regional  preferrence 
determined  by  magnetic  latitude.  Since  the  earth's  magnetic  field 
approximates  a  magnetic  dipole,  corrected  dipole  latitude  is  often  used 
instead  of  geographic  latitude.  The  auroral  zone  is  generally  accepted  to  be 
the  region  between  60  and  75  degrees  IN  Lat  (dipole  latitude),  on  average. 
This  is  the  most  likely  region  to  contain  the  auroral  oval ,  the  area  of  auroral 
activity.  Ann  Arbor,  Michigan's  dipole  latitude  is  about  53  deg  IN,  it’s 
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geographic  latitude  is  about  42  N.  The  IN  latitude  puts  Arm  Arbor  just 
outside  the  auroral  zone  . 

Figure  1.1  shows  the  auroral  oval  at  various  levels  of  magnetic 
disturbance. 


(o)  (b)  (c) 

Figure  1.1  Auroral  oval  for  various  levels  of  magnetic  disturbance  index  Q. 
Cases  (a),  (b),  and  (c)  are  for  Q=0,3,  and  7  respectively.  The  mean 
auroral  oval  is  similar  to  (b).  (after  Starkov  and  Feldstein,  1968). 

The  reason  for  the  auroral  zone's  existense  can  be  seen  when 
viewing  an  illustration  of  the  earth's  magnetic  field.  Figure  1.2  shows  the 
earth’s  magnetic  field  and  how  charged  particles  migrate  in  the  vicinity  of 
earth,  once  trapped  by  it's  magnetic  field. 

As  seen  in  the  illustration,  charged  particles  spiral  around  the 
magnetic  flux  lines  themselves,  while  drifting  around  the  earth.  When  the 
particles  reach  a  turning  point,  as  illustrated  in  the  figure,  they  reflect  back 
to  the  other  direction.  Thus  the  trapped  particles  (presumed  to  be  electrons 
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here)  oscillate  back  and  forth  between  reflection  points  as  they  drift  around 
the  globe,  to  a  first  approximation. 


spiral i ng  about 
magnetic 
flux  line 

Figure  1.2  Motion  of  trapped  particles  in  a  dipole-like  field,  (adapted  from 

Whalen  et  al.,  1974). 

These  electrons  are  very  energetic-  they  contain  many  keV's  of 
energy.  Compared  with  the  13.6  eV  it  takes  to  ionize  an  N2  molecule,  it  can 
be  seen  that  this  is  huge  amount  of  energy,  and  the  energy  from  a  single 
electron  can  be  distributed  to  thousands  of  ambient  particles. 

The  source  of  these  energetic  electrons  is  of  course  the  sun.  Solar 
flares,  giant  explosions  which  release  fiercely  hot  particles  from  the  sun, 
occur  more  frequently  during  high  sunspot  activity,  which  follows  an  11- 
year  cycle.  Some  of  these  energetic  particles  are  intercepted  by  earth's 
magnetic  field  and  may  become  trapped  by  the  region,  and  a  portion  of  them 
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are  steered  to  the  near-earth  region,  becoming  a  source  of  energy  for  the 
aurora.  Such  an  event  occurred  on  November  8,  1991 

1.1.3  The  November  1991  Measurement 

The  year  of  1991  was  an  active  year  for  the  star  which  we  call  the  sun. 
During  our  sun's  well  known  11  year  cycle,  1991  was  one  of  the  peak  years 
of  solar  activity.  On  November  8,  1991  solar  activity  was  fairly  high,  with  a 
10.7  cm  radio  wave  flux  (F10.7)  index  reading  of  200.4,  and  geomagnetic 
activity  near  earth  was  at  severe  storm  levels,  with  a  planetary  magnetic 
(Ap)  index  daily  average  at  65  on  November  8,  and  99  on  November  9.  On 
the  evening  of  the  8th  auroral  displays  were  witnessed  unusually  far 
equatorward,  as  far  south  as  Lubbock,  Texas  in  North  America.  In  Ann 
Arbor,  Michigan,  the  University  of  Michigan's  Near  Infrared  Spectrometer 
(NIRS)  was  operating  and  measured  the  auroral  display  overhead  from 
approximately  6  PM  to  12  PM  local  time,  in  the  spectral  region  of  8400  to 

o 

8900  A.  Frames  of  the  aurora  were  taken  at,  and  integrated  over,  30  second 
intervals  during  that  time. 

Figure  1.3  shows  the  measured  spectra  as  a  function  of  time.  The 
most  obvious  features  in  the  graph  are  the  more  intense  peaks  which  are 
attributable  to  atomic  oxygen  and  atomic  nitrogen  emission  lines.  The 
aurora  was  brightest  near  2000  hrs  local  standard  time  (LST),  or  8  PM,  then 
after  2100  LST  the  intensity  dropped  off  remarkably  and  remained  at  the 
much  less  intense  nightglow  levels  for  the  rest  of  the  measurement  period. 
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Figure  1.3  Measurement  of  aurora  taken  over  Ann  Arbor,  Michiga 


12  Purpose  of  Research 


The  purpose  of  this  research  was  to  analyze  the  data  and  determine 
which  atomic  and  molecular  species  were  responsible  for  the  emissions  in 
the  spectral  region  measured,  then  determine  species  brightnesses  and 
molecular  band  temperatures.  This  would  be  accomplished  by  comparing 
the  calculated  spectra  with  the  measured  spectra,  and  minimizing  the 
differences  between  the  two,  using  a  nonlinear  least  squares  method.  This 
method  of  analysis,  previously  having  been  applied  only  to  the  more  basic 
nightglow  spectra,  allows  for  accurate  determination  of  temperatures  and 
brightnesses,  limited  only  by  instrument  resolution.  The  method  also 
converges  to  the  solution  quite  rapidly,  requiring  fewer  calculations  than  a 
grid  search  or  other  numerical  method. 

Once  the  analysis  of  the  spectrometer  data  is  comi  ’ete,  a  study  of  the 
effects  the  auroral  electrons  had  on  the  upper  mesosphere  and  lower 
thermosphere  on  that  evening,  was  conducted.  Some  of  these  effects  are 
time  variation  of  vibrational  band  temperature  and  brightness,  correlation 
of  brightness  and  temperature  variation,  comparison  of  N  I  (non-ionized 
atomic  nitrogen)  emission  line  brightness  ratios  versus  expected  values, 
and  vibrational  population  of  N2  bands. 

1.3  Instrument  Description 

The  instrument  used  is  a  Visible  and  Near-infrared  Spectrometer 
(hereafter  abbreviated  "NIRS"),  developed  at  the  University  of  Michigan, 
and  is  a  1/4  meter  Ebert-Fastie  spectrometer  (Yee  et  al.,  1991).  The 
instrument  was  designed  to  be  used  on  rocket  flights,  but  was  modified  for 
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ground  based  operations.  A  brief  description  of  the  NIRS  instrument  and 
its  application  follows. 

The  NIRS  optical  components  are  illustrated  in  Figure  1.4,  and  will 
be  described  in  this  section.  The  complete  cooling  system  and  control 
electronics,  however,  will  not  be  discussed  here. 


Figure  1.4  The  University  of  Michigan  Near-infrared  Spectrometer. 

Telescope  and  Baffle.  NIRS  observes  the  atmosphere  through  a  two- 
stage  baffled  f/5  telescope.  The  telescope  provides  the  light  collection  and 
pointing,  and  directly  feeds  the  light  collected  to  the  spectrometer.  The 
baffle  is  similar  in  design  to  the  ones  used  in  the  Visible  Airglow 
Experiment  (VAE)  photometers  on  the  Atmospheric  Explorers  Satellites, 
and  the  Fabry-Perot  Interferometer  (FPI)  on  the  Dynamic  Explorer 
Satellites.  The  purpose  of  the  baffle  system  is  to  eliminate  stray  light  from 
entering  the  optics.  Internal  baffles  are  also  installed  to  further  minimize 
scattered  light. 
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Optics.  Optical  parameters  are  listed  in  Table  1.1.  The  optical  design 
was  fully  raytraced  and  carefully  optimized  to  minimize  the  off-axis 
blurring  due  to  optical  aberration. 


PARAMETER 

SPECIFICATION 

Spectral  range 

8355  -  8900  A 

Spectral  channels 

385  ( 1  used  for  dark  count) 

Wavelength  per  channel 

1.3983  A 

Spectral  resolution  (FWHH) 

12.6  A 

Telescope  f  -number 

f/5 

Objective  lens  focal  length 

150  mm 

Field  of  view 

1717  x  2.52  deg 

Input  slit  size 

0.1877x6.62  mm 

Grating  density 

600  lines/mm 

Interference  order 

1 

Detector  type 

Thompson  TH7882  CCD 

Detector  array  size 

388  x  288  pixels 

Detector  pixel  size 

23  x  23  mm 

Detector  temperature 

1 

o 

1 

*0 

n 

Sensitivity 

1.6  el/sec/R/A  @8600  A 

Table  1.1  NIRS  instrument  specifications. 

Ebert-Fastie  Spectrometer.  NIRS  uses  a  1/4  meter  Ebert  mirror  and  a 
600  lines/mm  holographic  grating,  and  operates  on  the  first  interference 
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order.  A  blocking  filter  which  eliminates  photons  with  wavelength  less 
than  6000A  is  installed  to  prevent  photons  with  higher  order  dispersion 
from  reaching  the  detector. 

Detector.  The  grating-dispersed  spectrum  is  imaged  onto  a  Charged- 
Coupled  Device  (CCD)  detector  using  on-chip  binning  of  photons.  The  CCD 
used  is  a  MOS  (metal-oxide  semiconductor)  optical  detector  composed  of 
111,744  (388  x  288)  independent  sites,  or  pixels,  where  photon-induced 
charge  is  stored.  The  number  of  charge  induced  is  porportional  to  the  light 
intensity  and  can  accumulate  over  an  extended  period  of  time  with  its 
exposure  length  (integration  time)  controlled  by  the  frequency  of  transfer 
between  the  storage  frame  and  image  frame.  The  integration  time  used  for 
this  study  was  30  seconds. 

The  CCD  detector  system  used  is  a  modified  Photometric  CCD 
camera  head  with  a  Thompson  TH7882  frame  transfer  CCD  chip  installed 
within  it.  The  chip  is  cooled  using  a  three  stage  thermo-electric  cooler  with 
flowing  liquid  coolant  in  order  to  decrease  the  thermal  noise.  This  thermal 
noise  exhibits  itself  in  a  CCD  as  dark  counts  (charge  induced  with  no  light 
entering  the  detector).  Usually,  a  CCD  temperature  increase  of  a  mere  7 
degrees  Celsius  can  double  the  thermal  noise  or  dark  counts.  Thus  in 
order  to  keep  thermal  noise  well  below  signal  levels  it  is  important  to  cool 
the  chip  as  much  as  possible.  The  temperature  was  predominantly  within 
-65  to  -75  degrees  Celsius  for  this  study. 

Sensitivity.  The  absolute  sensitivity  obtained  is  based  on  an  NBS 
(National  Bureau  of  Standards)  lamp  of  known  spectral  shape  and 
brightness.  Figure  1.5  gives  the  results  of  the  calibration  with  1  x  288  pixel 
CCD  binning.  As  can  be  seen  in  the  graph,  there  is  somewhat  less 
sensitivity  at  the  end  pixels  or  channels  due  to  the  grating  slightly 
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overlapping  the  CCD  chip.  Approximately  0.10  digital  counts/sec/R/A  tor 
1.6  el/sec/R/A)  sensitivity  is  obtained  in  the  spectral  region  near  86OOA  (the 
analog/digital  conversion  factor  is  15.6  electrons/digital  count). 


NIRS  ABSOLUTE  SENSITIVITY 


Wavelength  (A) 

Figure  1.5  Spectrometer  sensitivity  for  a  1  x  288  pixel  CCD  binning  pattern 
calibrated  January,  1991. 


Resolution.  The  NIRS  resolution  can  be  determined  from  the 
number  of  pixels  or  channels  it  takes  to  reduce  the  transmission  to  1/e  of 
the  reference  pixel's  transmission.  With  the  use  of  Figure  1.6,  it  can  be 
seen  that  1/e  from  the  peak  of  the  instrument  function  is  about  0.04  on  the  y- 
axis,  and  the  number  of  pixels  between  the  two  points  where  the  curve 
intersects  0.04  is  about  9.7  Thus  the  1/e  resolution  is  13.6  A. 

Another  measure  of  resolution  used  by  many  researchers  is  the 
measure  of  full  width  at  half  height  (FWHH)  of  ne  instrument  function. 
This  can  also  be  obtained  from  Figure  1.6.  The  full  width  at  half  height  for 
the  NIRS  instrument  is  12.6  A. 
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NIRS  INSTRUMENT  FUNCTION 


Figure  1.6  Instrument  function  for  the  spectrometer.  Shows  how  incident 
radiation  corresponding  to  channel  14  is  transmitted  to  surrounding 
channels,  thus  affecting  instrument  resolution. 


Chapter  2  will  show  how  data  from  the  instrument  was  analyzed  and 
how  the  atomic  and  molecular  parameters  were  recovered.  Chapter  3 
reveals  results  of  the  analysis  for  the  atomic  species  emissions,  and 
Chapter  4  gives  results  for  the  molecular  species.  Chapter  5  discusses 
conclusions  drawn  from  this  study  and  gives  recommendations  for  further 
studies. 
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CHAPTER  2 


ANALYTICAL  METHODS 


2.1  Non-linear  Least  Squares  Method 

The  nonlinear  least  squares  method  is  an  iterative  technique  similar 
to  Newton's  method.  It  requires  a  first  guess,  though  not  necessarily  close 
to  the  actual  solution,  then  in  each  iteration  it  finds  the  fastest  descent  to 
the  minimum  of  the  function  being  resolved.  A  similar  technique  was 
demonstrated  by  Killeen  and  Hays  (1984)  for  the  analysis  of  Fabry- Perot 
interferometer  data. 

A  remarkably  few  number  of  iterations  is  required  for  convergence  to 
occur,  and  the  criterion  for  convergence  can  be  arbitrarily  chosen  by  the 
user.  In  this  study  a  value  of  0.5%  of  each  parameter  was  used  as  the 
convergence  criterion,  a  fairly  strict  requirement.  The  number  of  iterations 
can  be  as  few  as  3  or  4  depending  on  the  convergence  criterion  and  goodness 
of  lit.  A  description  of  the  least  squares  method  used  follows. 

As  an  introduction  to  notation  used  in  the  analysis,  the 
measurement  will  be  designated  as  Y  ,  and  f(x)  is  a  known  function  or  set 
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of  functions,  called  a  fitting  function,  which  we  can  use  to  predict  the  value 
of  Y  by  varying  x  in  an  appropriate  manner. 

To  fit  more  than  a  single  parameter  requires  a  greater  number  of 
independent  simultaneous  measurements,  or  channels,  than  fitting 
parameters,  and  the  greater  the  number  of  channels  over  the  number  of 
fitting  parameters  the  better.  Thus  the  number  of  degrees  of  freedom  is  the 
difference,  N-n,  where  N  is  the  number  of  instrument  channels,  and  n  is 
the  number  of  parameters  used  in  the  fit. 

ty 

The  objective  is  to  find  the  fitting  parameter  x,  by  minimizing  %~- 
This  is  done  by  finding  the  weighted  difference  between  measurement  and 
prediction: 


r  = 


N-, 


■X -Hiw.wf 


*=i 


2.1 


where  crk  is  the  measurement  error  of  the  kth  channel,  and  N  and  n  are  as 
defined  above. 

2.1.1  Solving  for  Parameters 

Performing  a  first  order  Taylor  expansion  of  equation  2.1,  about  a 
point  a,  of  the  function  f,  we  get: 


r  = 


N-n 


L 


<>) 


2.2 


As  seen  in  equation  2.2,  if  the  instrument  has  more  than  one  data  channel, 
then  the  measurement  output  Y,  the  output  of  the  fitting  function  f,  and  the 
derivatives  of  f,  all  must  be  vectors,  with  the  number  of  dimensions  of  the 
vectors  being  equal  to  the  number  of  instrument  channels.  This  also 
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requires  some  manipulation  of  the  output  of  the  fitting  function  (through 
the  use  of  a  convolver),  so  that  a  comparison  of  the  fitted,  synthetic 
spectrum,  with  the  measurement,  on  a  channel-by-channel  basis  can  be 
accomplished.  More  on  this  will  be  discussed  in  section  2.3. 

If  the  fitting  function  requires  more  than  one  variable,  we  need  to 
write  equation  2.2  in  multi-variable  form,  so  x  and  a  become  multi¬ 
dimensional,  or  vectors,  and  are  labelled  as  such.  Equation  2.2  now 
becomes: 


r  = 


_L_v_L  ly 

N  -  n  ^  <7k  j  k 


-/*(")“ 


v  -LSEI 
k  9x, 


2.3 


where  i  varies  from  1  to  the  number  of  parameters  being  analyzed,  n,  and  k 
varies  from  1  to  the  number  of  channels,  N. 

We  can  find  the  best  possible  value  of  each  variable,  x,,  which  can 

represent  species  brightness  or  temperature  for  example,  by  minimizing 
the  above  expression.  To  minimize  X'  we  require  the  relations: 


2.4 


To  simplify  matters  let's  assume  we  have  only  two  channels  (N= 2) 
and  two  parameters  to  fit  (n  =  2).  Since  we  need  more  charm-ds  than 
parameters,  we  know  this  would  not  work  when  trying  to  calculate  x~ 
(because  of  a  divide  by  zero  error),  but  the  system  of  equations  can  still  be 
solved,  so  let's  continue  anyway  to  illustrate  the  mathematics.  Now 
equation  2.3  becomes: 
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2.5 


i/  ^  dfx(a) 


U,  - a ,) 


1  L  x  N  #■»(«),. 

+ —  <  K,  -/,(«) - 5 - ( Xj  -a,) - - - (*2  -«,) 

a,  *  eh,  r/r. 


Using  equation  2.4,  neglecting  third  order  terms,  and  abbreviating  /,  2(ct)  as 
/[  2,  it  can  be  shown  that: 


cot,  ctx,  cr,  ^  ctx,  j  a2  dx2 


and: 


=  -2~-  My,  -a,)]-  —  =  <>■ 

eCc2  r/x2  <72  \  "  '  »X2  "  y  (7,  e/x, 


Simplifying  and  rearranging  equation  2.6  becomes: 


1  f  ^  .  .  1  df\  df\  1  dfx  r  , 

-  a, )  +  —  -  - (x2  -^2)  =  — —  y,  -  /, (a ) 

<7,  ^  ox,  )  <J2  ox,  dx2  <7,  ox, 

1  df2  df2  ,  I  ( <?/,  V,  .  ..  ,  1  df2  {v  r 

(X|  f/| )  +  ( Xi  a-,)—  i y,  ./\(^)J 

cr,  ox,  ox,  cr,  ^  ox2  y  *  ‘  cr2  dx2 


which  in  two-dimensional  matrix  form  looks  like: 


ffi  )  _  . 

v  r?x,  J  c?x,  f)x2  ct,  1  ^ 

i^Y  — u2-«2) 

<9x,  r)x2  \  flx2  )  L  - 


-vL(i,i -/,(«)) 


(7,  r?X, 


(7n  r?x, 
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Thus  in  general  we  get  n  equations  and  n  unknowns,  and  we  are  free 
to  choose  n,  the  number  of  fitting  parameters,  as  large  as  we  wish,  up  to  the 
number  of  channels  on  the  spectrometer  (the  number  of  degrees  of  freedom) 
minus  one.  The  matrices  in  2.9  are  of  the  form  Ax  -  B  and  we  wish  to  solve 
for  the  x  matrix.  After  doing  so  we  can  add  the  a-  s,  which  are  the  initial 
values,  to  arrive  at  the  x's.  One  can  use  an  LU  decomposition  or  other 
technique  to  solve  for  the  x  matrix,  but  here  we  will  use  the  simple  matrix 
inversion  method:  x  =  A~'B  (code  adapted  from  Bevington,  1969),  since  the 
inverve  matrix  elements  will  come  in  use  later  on  when  determining  the 
most  probable  errors. 

This  method  very  easily  becomes  an  iterative  method  by  reentering 
the  loop,  substituting  the  u/s  with  the  .r/s  obtained  from  the  previous 
iteration,  or  if  it's  the  first  iteration  arbitrarily  set  the  o'  s  to  some  initial 
guess.  In  one  line  form  the  iteration  formula  may  be  written: 

k 

where  Xj  mis  the  i th  parameter  from  the  mth  iteration,  A",.1  is  the  inverse  of 

matrix  A,  and  is  also  a  square  matrix,  and  B  has  N  rows  and  1  column. 
We  are  now  ready  to  test  to  see  when  the  iteration  can  be  stopped. 

2.1.2  Iteration  Testing 

Since  we  have  the  iteration  formulae,  we  now  need  to  know:  1)  when 
to  stop  the  program  from  looping,  and  2)  how  good  the  calculated  spectra 
fits  the  measured  spectra. 

The  first  question  can  be  answered  by  simply  looking  at  the  relative 
change  of  each  parameter  after  each  iteration.  If  the  percentage  change  of 
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each  parameter  is  less  than  some  arbitrarily  small  number,  then  a  simple 
"IF"  statement  in  FORTRAN  can  stop  the  loop.  In  this  study,  if  the  relative 
change  of  each  parameter  was  less  than  0.5%  then  the  loop  was  halted,  and 
convergence  has  occurred.  An  exception  to  this  rule  was  made  for  those 
parameters  whose  intensity  was  less  than  100  R,  which  generally  indicated 
a  poor  fitting  parameter  to  that  particular  spectra. 

The  question  of  how  good  is  the  fit  can  be  answered  in  part  by  looking 
at  the  value  of  %2  ■  For  a  perfect  fit,  x~  =1*  If  the  value  is  much  greater  than 
say  5,  then  the  synthetic  spectrum  does  not  fit  the  measured  spectrum  well 
enough  to  say  the  spectrum  is  fully  analyzed. 

Now  that  we  have  the  method  to  optimize  any  variable  we  choose,  we 
need  to  know  what  is  in  the  function  f  ,  which  we  will  use  in  equation  2.9  to 
determine  the  A  matrix  which  will  then  be  inverted.  The  function  f 
inputs  all  parameters  used  in  the  analysis  and  outputs  a  channel  by 
channel  intensity  distribution  for  use  in  comparing  with  the  instrument 
data. 

To  do  this  requires  two  basic  steps:  1)  take  input  parameters  (species 
brightnesses  and  temperatures)  and  determine  a  theoretical  emission 
spectrum,  2)  take  emission  spectra  and  calculate  the  signal  received  for 
each  wavelength  interval  corresponding  to  the  NIRS  instrument  channels, 
and  sum  the  individual  intensities  to  get  a  total  intensity  for  each  channel. 
Step  1  requires  a  band  model  for  each  molecular  transition  being 
investigated,  and  step  2  requires  a  convolver  which  depends  on  individual 
instrument  characteristics. 
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2.2  Band  Models 


As  mentioned  before,  a  band  model  is  needed  for  each  molecular 
transition  being  investigated.  Due  to  the  rich  variety  of  emissions  in  the 
auroral  measurements,  this  study  required  several  band  models.  The 
molecular  bands  used  to  fit  the  aurora  data  are  listed  in  Table  2.1. 


Band  Name 

Electronic 

Transition 

Vibrational 

Transitions 

(v'.v") 

Model  Author 

N2  1st  Positive 

(1,0), (2,1), (3,2), 
(4,3) 

R.L.  Gattinger 

O2  Atmosperic 

(0,1), (1,2), (2, 3) 

P.B.  Hays 

OH  Meinel 

N/A 

(6, 2), (7, 3) 

J.H.  Yee  (1988) 

N2  Infrared 
Afterglow 

B"ru-Bynx 

(6,2),(11,6) 

(12,7) 

E.J.  Bucsela, 
T.J.  Smith  (1992) 

N*  Meinel 

A2nm-x% 

(4, 2), (5, 3) 

E.J.  Bucsela 
(1991) 

Table  2.1  Molecular  bands  and  band  models  used  to  analyzed  the 

November  8,  1991  aurora  measurements.  The  OH  Meinel  bands  are 
due  to  pure  vibration-rotation  transitions,  so  an  electronic  transition 
is  not  applicable  for  this  case. 


The  basic  principle  of  each  of  the  above  models  is  the  same:  input  the 
desired  vibrational  transition,  temperature,  and  band  brightness,  and  the 
model  outputs  an  array  of  emission  line  wavelengths  and  an  array  of 
associated  line  strengths.  The  number  of  lines  output  is  on  the  order  of 
hundreds  to  over  a  thousand  for  each  band.  An  example  of  output  from  a 
band  model  is  given  in  Figure  2.1.  Figure  2.2,  shown  next  to  Figure  2.1  for 
comparison  purposes,  is  what  Figure  2.1  looks  like  once  the  lines  are 
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convolved,  i.e.,  binned  by  wavelength  and  converted  to  a  "measurement"  as 
seen  by  the  spectrometer.  The  convolver  subroutine  is  discussed  in  the  next 
section. 

Once  all  the  appropriate  models  have  been  called  in  the  computer 
program,  the  arrays  of  wavelengths  and  strengths  are  stacked  so  that  you 
are  left  with  one  array  of  wavelengths  and  one  array  of  line  strengths. 
These  arrays  are  next  fed  into  the  convolver  subroutine. 

2.3  NERS  Convolver 

The  NIRS  Convolver  subroutine  was  written  by  J.H.  Yee  for  use  in 
determining  atomic  emission  intensities  and  molecular  emission 
temperatures  and  intensities  from  the  instrument  data.  I  will  merely 
describe  the  principles  used  to  write  the  subroutine. 

The  idea  of  the  convolver  subroutine  is  simple:  take  the  arrays  of  line 
wavelengths  and  associated  line  strengths,  compute  the  signals  according 
to  the  instrument  function  and  sensitivity  profiles  shown  in  Figures  1.3  and 
1.4,  then  add  up  the  signals  within  wavelength  bins  which  correspond 
to  instrument  channels.  This  procedure  takes  thousands  of  individual 
lines  and  converts  them  to  a  384  channel  predicted  measurement  which 
can  then  be  used  to  compare  with  the  actual  data.  An  example  of  a 
convolved  molecular  band  was  shown  in  fig.  2.2. 

It  should  now  be  evident  that  the  function  f  sought  in  section  2.1  is 
simply  a  summation  of  all  lines,  molecular  band  lines  and  atomic  lines, 
convolved  into  a  384  channel  prediction.  Thus  we  can  build  a  subroutine 
which  calls  the  desired  band  model  subroutines  with  input  band 
temperatures  and  intensities,  then  convolves  the  many  lines  into  a  useful 
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Brightness  (counts/30  sec) 


8500  8550  8«X)  8650  o87(X)  8750  88(X) 

Wavelength  (A) 

Figure  2.1  An  example  of  band  model  output  -  N2  first  positive  (2,1)  band  at 
1000  R  band  intensity  and  300  K. 


8500  8550  8600  8650  §700  8750  8800 

Wavelength  (A) 

Figure  2.2  Convolved  N2  1st  positive  (2,1)  band  at  1000  R  band  intensity  and 
300  K. 
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form.  This  subroutine  is  like  the  function  f  ,  used  earlier  in  this  chapter, 

which  we  can  use  to  find  the  f(a)  or  —  arrays  in  equation  2.9. 

dxt 


2.4  Error  Analysis 

When  measuring  the  aurora  with  a  CCD  spectrometer,  there  are 
three  basic  types  of  noise  which  must  be  dealt  with  in  some  way:  1)  thermal 
noise  (dark  counts),  2)  instrument  bias  (read  noise),  and  3)  photon 
statistics.  The  treatment  of  each  of  these  will  be  discussed  separately. 

2.4.1  Thermal  Noise 

A  CCD  detector  produces  thermal  noise  or  dark  counts  which  is 
greatly  dependent  on  temperature.  As  mentioned  previously  the  detector 
system  used  for  this  study  incorporated  a  385  channel  CCD  binning  pattern, 
with  384  columns  of  pixels  used  as  data  channels,  and  4  columns  of  pixels 
used  to  measure  dark  counts. 

The  dark  count  columns  of  pixels  were  covered  so  that  no  external 
radiation  would  be  detected  at  these  columns  of  sites,  which  were  used  as  a 
single  channel.  It  is  assumed  that  any  counts  detected  from  this  channel 
must  be  due  to  thermal  noise  which  is  dependent  on  the  CCD  temperature, 
assumed  to  be  uniform  across  the  chip.  Using  these  assumptions  the 
thermal  noise  can  be  calculated  and  carefully  removed,  as  much  as 
possible,  by  simply  subtracting  the  number  of  dark  counts  detected  by  the 
385th  channel  from  each  of  the  other  384  channels.  One  way  to  do  this  is  by 
using  an  approximation  for  the  CCD  temperature  from  the  dark  count 
measurement: 
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2.11 


Tnv{°C)  =  -50  -  //tRJn(-^-  AO, 

where  H  and  Cl  are  simply  coefficients  for  the  385th  channel  (the  dark 
channel)  as  measured  during  calibration,  C  is  the  number  of  dark  counts 
measured  during  the  30  second  integration  period  (At),  and  -50  degrees  was 
chosen  arbitrarily  as  a  reference  temperature.  To  subtract  the  number  of 
dark  counts  from  the  i th  channel  measurement,  we  use  the  relation: 

c,  =  -  Ar  •  Cl, cxp(^O) ,  2.12 

where  a<>  is  the  original  measurement  of  the  i th  channel,  and  c<  is  now 
the  measurement  adjusted  for  thermal  noise.  It  is  agreed  there  are  a  few' 
redundant  calculations  in  the  above  equations,  but  these  equations  were 
chosen  for  ease  of  understanding,  versus  ease  for  computation. 

2.4.2  Instrument  Bias  (Read  Noise) 

This  noise  is  introduced  by  the  readout  amplifier  of  the  CCD 
electronics.  For  each  read  that  occurs,  a  slight  error  in  the  number  of 
counts  read  is  introduced.  The  magnitude  of  the  error  introduced  is  on  the 
order  of  1  count/read,  where  a  read  must  be  accomplished  to  retrieve  data 
from  each  channel,  for  each  frame  of  data.  The  instrument  bias  is 
corrected  for  by  simply  taking  the  bias  (determined  from  laboratory 
calibration)  and  subtracting  it  from  the  measurement  of  each  channel. 

2.4.3  Photon  Statistics 

This  type  of  noise  cannot  be  corrected  for  as  with  the  above  types  of 
errors  because  we  are  uncertain  of  the  amount  of  the  error  involved  in  each 


25 


measurement  within  known  bounds.  However,  statistical  errors  usually 
follow  a  Poisson  distribution: 


where  I  is  the  measured  signal,  expressed  in  units  of  number  of  electrons. 

From  this  relation  we  can  see  that  the  greater  the  signal,  the  smaller 
the  relative  error.  To  exploit  this  all  we  need  to  do  is  increase  the 
instrument’s  integration  time,  or  time  the  shutter  remains  open,  to 
increase  intensity,  thereby  lowering  relative  error. 

As  an  example  of  how  this  integration  helps,  let's  compare  the 
relative  error  of  low  and  high  intensity  measurements  at  1  versus  30  second 
integration  time.  During  the  early  evening  hours  (when  the  aurora  was 
active)  the  spectrometer  recorded  counts  between  10  and  100  counts  per 
second.  Once  the  aurora  quieted  down,  after  9  PM,  nightglow  emissions 
predominated  and  only  about  1-10  counts  per  second  were  recorded. 

Table  2.2  shows  a  comparison  of  the  statistical  errors  calculated  from 
equation  2.13,  remembering  the  error  is  based  on  electron  counts  not  digital 
counts.  Thus,  the  /  in  equation  2.13  is  the  number  of  digital  counts 
measured,  multiplied  by  15.6  electrons/digital  counts,  which  is  then 
substituted  into  the  equation  for  error  determination.  As  an  added  note, 
throughout  the  remainder  of  this  writing,  "counts"  will  refer  to  digital 
counts,  not  number  of  electrons.  Conversion  to  number  of  electrons  is 
necessary  when  figuring  photon  statistics. 

As  can  be  seen  from  the  table,  a  30  second  integration  period 
decreases  statistical  error  by  a  factor  of  about  5!  Of  course  this  increase  in 
accuracy  is  not  free.  Since  we  are  integrating  over  30  seconds  we  are  losing 


that  much  temporal  resolution,  but  not  enough  to  be  deemed  important.  A 
30  second  resolution  is  sufficient  for  most  any  study. 


Also  seen  from  Table  2.2  is  the  fact  that  measurements  taken  during 
the  auroral  maximum  period,  30  second  integration  time,  are  accurate  to 
within  0.5%  (for  peak  intensity  channels)  to  2.5%  (for  least  intense 
channels).  However,  between  9-12  PM,  measurements  are  accurate  only  to 
within  1.5%  to  8.0%. 


Measurement 

Condition 

Intensity 

Relative  Error 

(counts/sec) 

1  sec 
int  per 

30  sec 
int  per 

Low  int. 
nightglow 

1 

25% 

4.6% 

High  int.  nightglow. 
Low  int.  aurora 

10 

8.0% 

1.5% 

High  int.  aurora 

100 

2.5% 

0.5%, 

Table  2.2  Comparison  of  statistical  error  ol  measurements  with  1  second 
versus  30  second  integration  period  under  different  brightness 
conditions. 


An  estimate  of  statistical  error  of  a  particular  fitting  parameter  can 
be  obtained  quite  simply.  The  uncertainty  in  fitting  they  th  parameter, 

is  given  by: 


2.14 


where  Ejj  is  the  j  th  diagonal  element  of  the  error  ma  ’rix.  r.  The  error 
matrix  is  basically  the  inverse  matrix,  A'1,  for  which  we  have  already 
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solved  to  recover  the  a/s.  The  a's,  then,  are  the  analytical  errors  of  the 
fitting  parameters. 


2.4.4  Systematic  Error 

The  predominant  type  of  systematic  error  in  the  spectrometer 
measurements  occurs  because  of  absolute  sensitivity.  The  sensitivity 
calibrations  are  accomplished  using  a  known  continuum  source,  with  its 
brightness  traceable  to  an  NBS  lamp,  as  described  in  chapter  1  Calibrating 
an  instrument  to  this  lamp  does  not  guarantee  absolute  accuracy.  In  fact, 
the  NIRS  instrument  can  only  claim  an  accuracy  of  15%  due  to  the 
uncertainty  of  the  standard. 

Other  types  of  systematic  error  arise  because  of  imprecise  methods 
used  in  various  stages  of  the  fitting  function.  Some  examples  are 
assumption  of  only  electric  dipole  transitions  in  band  models,  imprecise 
molecular  constants  used  to  calculate  band  model  output,  assumption  of 
thermal  equilibrium  among  rotational  levels,  etc. 

2.5  Preliminary  Results 

An  analysis  of  the  spectrometer  data  was  done  in  the  manner 
described  in  previous  sections.  Table  2.3  lists  the  parameters  used  in  the 
non-linear  least  squares  technique  to  create  a  synthetic  best-fit  spectrum. 
The  table  of  parameters  is  just  one  of  nearly  600  frames  analyzed  and 
represents  typical  intensities  and  temperatures  recovered  during  the  period 
of  higher  auroral  activity. 
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On  page  31,  Figure  2.3  shows  how  those  parameters  fit  the  spectrum. 
A  two-parameter  continuum  was  used  to  allow  for  better  fitting  of  the  data 
since  it  is  presumed  that  the  background  radiation  is  not  a  constant  across 
the  spectrum.  At  least  10  atomic  oxygen  and  nitrogen  lines  were  found  to 
exist  in  the  spectrum  and  14  molecular  bands  were  analyzed  as  well.  A 
detailed  description  of  the  analysis  of  atomic  species  is  given  in  chapter  3, 
and  the  description  of  the  molecular  species  analysis  is  given  in  chapter  4. 
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Value 

Units 

Error* 

8.51 

R 

1.69 

25  87 

R 

0.63 

8548.51 

R 

18.05 

!ME£59i 

R 

8.01 

434.60 

R 

11.10 

1131.00 

R 

24.05 

521.50 

R 

19.79 

4869.28 

R 

34.42 

5627.37 

R 

35.83 

804.16 

R 

37.43 

1164.04 

R 

23.67 

707.81 

R 

33.89 

32834.50 

R 

273.45 

24009.21 

R 

334.62 

3401.74 

R 

80.47 

490.62 

R 

500.01 

383.05 

K 

4.08 

11115.95 

R 

97.33 

432.07 

K 

3.17 

4559.97 

R 

162.90 

1356.45 

R 

64.53 

637.69 

K 

16.17 

1703.79 

R 

58.68 

2501.85 

R 

mnumm 

158.84 

K 

3.70 

4103.61 

R 

150.08 

1438.56 

R 

114.62 

5919.23 

R 

439.84 

640.40 

K 

30.44 

3186.54 

R 

141.76 

1589.56 

R 

196.00 

Parameter 


Hi  1  continuum 


Lo  1  continuum 


0  1  8446.5  A 


N  I  8567.7  A 


(12,7) 


P 


4,2 


Time:  19.0167 


=  21.67 


Table  2.3  List  of  the  31  parameters  used  in  this  study  to  analyze  a  typical 
measurement  spectrum,  with  an  integration  period  of  30  seconds. 

(*  Note  that  since  the  value  of  Chi-squared  is  not  equal  to  1  the  error 
values  in  the  above  table  are  not  representative  of  the  true  analytical 
uncertainty). 


AURORAL  EMISSION  COMPOSITION 


oas  oe/xN3)  axisn3Xni 


31 


Figure  2.3  Total  analyzed  composition  of  the  November  8,  1991  aurora  over 
Ann  Arbor,  Ml.  The  top  solid  line  indicates  the  measurement,  dotted 
indicates  fitted  spectrum,  lower  lines  indicate  fitted  molecular 
parameters  and  continua. 


CHAPTER  3 


ATOMIC  EMISSIONS 


3.1  Atomic  Structure 

Before  the  features  of  atomic  emissions  can  be  discussed,  atomic 
structure  must  first  be  understood.  I  will  not  go  into  great  detail  of  the 
structure,  but  enough  to  provide  a  background  so  that  an  understanding  of 
the  discussion  that  follows  is  possible.  This  section,  as  well  as  a  much 
more  detailed  discussion  can  be  found  in  Herzberg  (1945). 

The  Schroedinger  equation,  which  describes  the  motion  of  the 
electrons  in  an  atom,  requires  at  least  four  quantum  numbers  to  uniquely 
define  the  state  of  a  particular  atomic  electron: 

n:  principle  quantum  number 
1:  orbital  quantum  number 
m,:  orbital  magnetic  quantum  number 
ms:  spin  magnetic  quantum  number. 

According  to  the  Pauli  Exclusion  Principle,  no  two  electrons  can  have  the 
same  set  of  four  quantum  numbers  in  an  atom  or  molecule. 


The  principle  quantum  number,  n,  is  related  to  the  probability  of 
finding  the  electron  a  certain  distance  from  the  nucleus.  For  historical 
reasons,  all  electons  with  the  same  principle  quantum  number  are  said  to 
be  in  the  same  electonic  shell,  identified  by  the  letters:  K,L,M,N,0,  or  P  for 
n=l,2,3,4,5,  or  6,  respectively.  Likewise  electrons  having  the  same  value  of 
both  n  and  1  quantum  numbers  are  said  to  be  in  the  same  subshell, 
designated  as:  s,p,d,f,g,  or  h  when  1=0, 1,2, 3, 4,  or  5  respectively. 

As  an  example  of  notation,  the  ground  (or  valence)  electronic 
configuration  of  atomic  oxygen  is:  Is2  2s7  2  p4.  This  means  there  are  2 
electrons  in  the  Is  subshell  (or  K  shell,  n=l).  Likewise  the  2s  subshell  has  2 
electrons,  and  2p  has  4  electrons.  The  four  quantum  numbers  are  related 
as  given  in  Table  3.1. 


Quantum  Number 

#  Possible  Values 

Range  of  Values 

n 

1 

1 

n 

mi 

21+1 

-l.-l+l . 1 

ms 

2 

-1/2,  1/2 

Table  3.1  Relationship  of  the  four  atomic  quantum  numbers. 

The  presence  of  an  external  magnetic  field  gives  rise  to  splitting  in 
spectral  lines,  because  of  the  field's  interaction  with  the  spin  of  the  electron 
(anomalous  Zeeman  Effect).  This  splitting  can  be  very  small  or  quite 
measureable.  An  example  of  this  can  be  seen  in  Table  3.2.  The  splitting  of 
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the  three  01  lines  is  on  the  order  of  0.1A,  but  for  the  NI  8  series,  the 
splitting  is  on  the  order  of  10's  of  Angstroms. 

To  arrive  at  energy  states  of  an  atom  we  must  add  the  respective 
angular  momenta  of  the  individual  electrons.  The  energy  state,  or  term 
type,  of  an  atom  is  written  in  the  general  form: 


where  L  is  the  resultant  angular  momentum  quantum  number,  S  is  the 
resultant  spin  quantum  number,  and  J  is  the  resultant  total  angular 
momentum  quantum  number,  and  can  take  on  values  from 

J  =  |L  +  S| . |Z--5|.  When  L=0,l,2,  or  3  for  an  atom,  the  state  is  denoted  as 

a/an  S,P,D,  or  F  state  (analagous  to  the  electron  shell  designation). 

The  ground  state  of  atomic  oxygen,  l.v22.r2//,  has  term  type: 


From  this  notation  we  know  that  L=1  (P),  S=1  (2S+1=3),  and  since  in 
general,  a  state  with  a  left  superscript  of  three  denotes  a  triplet,  or  energy 
state  split  into  three  states,  J  can  take  on  a  value  of  either  2,  1,  or  0.  The 
exception  to  the  rule,  that  the  superscript  equals  the  number  of  split  levels, 
occurs  when  the  resultant  quantum  number,  L,  is  less  than  the  resultant 
quantum  number  S.  In  that  case  the  number  of  split  levels  is  equal  to  L. 
We  can  now  go  into  specifics  of  atomic  oxygen  and  nitrogen  emissions  and 
what  causes  them. 
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3.2  Atomic  Lines  in  the  Aurora 


There  are  a  number  of  atomic  lines  in  the  8400  8900A  spectral  region 
which  must  be  included  in  any  analysis  of  auroral  emissions.  The  atoms 
that  need  to  be  considered  are  primarily  those  of  the  O  and  N  atoms,  since 
these  atoms  are  some  of  the  most  abundant  in  the  altitude  region  of  the 
aurora,  100-200  km.  A  more  complete  analysis  would  need  to  include 
emissions  from  N*  as  well.  The  ionic  species  were  neglected  for  this  study 
because  it  was  attempted  to  keep  the  number  of  parameters  to  a  minimum, 
and  their  emission  intensities  were  expected  to  be  small. 

3.2.1  Production  Mechanisms 

The  production  of  the  excited  states  of  the  O  and  N  atoms  responsible 
for  the  emissions  may  be  due  to  inelastic  collisions  of  the  form, 

e  +  X  — >  X  *  +e,  3.3 

or  the  result  of  atom-atom  or  molecule-atom  energy  transfer, 

M*  +  X  ->  M  +  X*,  3.4 

where  X  may  be  either  an  N  or  O  atom,  M  may  be  any  atom  or  molecule, 
and  an  asterisk  denotes  an  excited  species  (Jursa,  1985).  Another  possibly 
equal  production  mechanism  for  excited  N,  suggested  by  Fillipini,  et  al. 
(1982),  is: 


N,  +e-^N:‘  +  2e 
N:*  ->  N*  +N 
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Once  the  excited  species  of  0  or  N  is  created,  a  radiative  transition 
may  occur,  and  the  speed  of  transition,  or  the  inverse  of  radiative  lifetime, 
is  determined  by  selection  rules  of  the  transition.  For  a  detailed  discussion 
of  selection  rules  and  types  of  transitions,  again  see  Herzberg’s  book. 
Magnetic  dipole  and  electric  quadrupole  transitions  occur  considerably 
slower  than  electric  dipole  transitions.  Typical  radiative  lifetimes  for  the  N 
species  electric  dipole  allowed  emissions  are  10  seconds-  very  fast. 

3.2.2  Emission  Lines 

A  list  of  ail  known  01  (neutral  0),  011  (O4),  NI,  and  Nil  atomic  lines 
in  the  8400-8900A  spectral  region  is  given  in  Table  3.2.  The  numbers  in 
italics  indicate  lines  used  in  this  study,  and  the  three  8446A  01  lines  were 
combined  as  discussed  in  the  next  section.  A  brief  discussion  of  each 
atom's  excited  state(s)  as  well  as  findings  in  the  auroral  analysis  arc  given 
in  sections  3.3  and  3.4. 

The  O  and  N  lines  selected  from  table  used  to  analyze  the  aurora 
spectra,  are  primarily  the  strongest  lines,  indicated  by  the  relative  intensity 
in  Table  3.2.  The  less  intense  lines  (those  less  than  500  relative  intensity) 
and  ionic  species  were  omitted  to  minimize  parameters  involved. 
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Wavelength 

(A) 

Origin 

Relative 

Intensity 

Einstein 

Probability 

A(10*sec  ') 

1 

8426.16 

OI 

120 

2 

8438.74 

Nil 

550 

3 

8446.25 

OI 

810 

4 

8446.36 

1000 

5 

8446.76 

" 

935 

6 

8567.74 

NI  8 

500 

0.0458 

mm 

8594.00 

ii 

570 

0.190 

8 

8629.24 

n 

650 

0.238 

9 

8655.89 

n 

500 

0.099 

10 

8676.08 

Nil 

220 

11 

8680.28 

NI  1 

700 

0.191 

12 

8683.40 

" 

650 

0.133 

13 

8686.15 

" 

500 

0.079 

14 

8687.43 

Nil 

110 

15 

8699.00 

" 

110 

16 

8703.25 

NI  1 

500 

0.171 

17 

8710.54 

Nil 

160 

18 

wmmsmm 

NI  1 

570 

19 

8718.83 

ii 

500 

0.054 

20 

8820.43 

OI 

325 

21 

8728.89 

NI  1 

250 

0.030 

22 

8747.36 

•• 

200 

0.0079 

Table  3.2  Known  oxygen  and  nitrogen  atomic  emission  lines  from  8350- 
8900A.  Numbers  in  italics  indicate  atomic  lines  used  in  this 
study.  (Wavelengths  and  relative  intensities  from  Weast,  1988; 
transition  probabilities  from  Wiese  et  al.,  1966). 


3.3  Atomic  Oxygen 


3.3.1  Source  of  Emission 

The  ground  (valence)  state  electronic  configuration  of  the  O  atom  is 
\sl2s12p* .  As  seen  from  Figure  3.1,  the  source  of  the  01  (neutral  oxygen) 

o 

8446. 5A  emission  line  is  the  transition  between  Rydberg  states 
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(\s12s22pi)3p  3P  to  3.v  ySn .  Since  the  upper  state  CP)  is  a  triplet  multiplet, 
and  the  lower  state  is  a  singlet  multiplet,  there  are  actually  three  emission 
lines  with  slightly  different  wavelengths  coming  from  this  transition: 
8446.25,  8446.36,  and  8446. 76A,  with  relative  strengths  of  810,1000,  and  935 
respectively  (Weast,  1988). 


3.3.2  Analysis 

Since  the  instrument  has  a  FWHH  resolution  of  12. 6A  (from  section 
1.3.1),  it  is  not  able  to  distinguish  any  of  the  three  01  8446  lines  discussed  in 
the  above  section  apart  from  one  another,  and  we  cannot  determine  the 
relative  intensities  of  the  three  individual  lines.  Instead,  a  single  line  at 

o 

8446. 5A  was  used,  and  the  intensity  of  this  line  was  obtained  from  the  least 
squares  technique  discussed  in  section  2.1. 


!o"  r!  "o  'v:  *** 

0*  2S*  2oVo°) 


Figure  3.1  Energy  level  diagram  of  the  0  atom.  Solid  lines  correspond  to 
lines  observed  in  aurora;  the  dashed  lines  indicate  doubtful  or 
predicted  emissions.  The  numbers  in  brackets  are  multiplet 
numbers  from  Moore  (1945).  (from  Vallance  Jones,  1974). 
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Figure  3.2  shows  how  the  intensity  of  the  OI  8446  line  varied  with 
time  from  18.5  up  to  21.3  EST,  when  the  aurora  all  but  ended  over  Ann 
Arbor,  Michigan. 


O  I  8446  INTENSITY 


Local  Time 


Figure  3.2  Intensity  of  the  OI  8446A  line  versus  time  (EST). 

As  can  be  seen  from  the  graph,  the  8446  line  intensity  varies  quite 
rapidly  and  appears  to  be  very  sensitive  to  energetic  particle  flux.  From 
figure  3.1  we  know  that  the  0  atom  needs  at  least  11  eV  of  energy  gained 
from  collisional  processes  of  equation  3.1  or  3.2.  From  equation  3.1,  the 
electrons  in  the  11  eV  energy  range  are  secondary  electrons,  which  were 
created  from  other  more  energetic  secondary  electrons  or  primary 
electrons.  One  possible  candidate  for  a  particle  of  the  type  in  equation  3,2  is 
energetic  0  or  G+  created  from  the  aurora’s  influx  of  energetic  electrons. 
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3.4  Nitrogen 


3.4.1  Source  of  Emission 

The  number  of  atomic  nitrogen  lines  is  much  greater  than  that  of 

o 

atomic  oxygen  in  the  range  of  8400  -  8900  A.  This  is  simply  due  to  energy 
states  of  the  N  atom  being  in  the  right  place  and  of  the  right  separation  to  be 
present  and  distinct  from  each  other  in  this  spectral  region.  Figure  3.3 
shows  the  two  transitions  responsible  for  the  NI  1  and  8  series  of  emission 
lines. 


DC.  Sit’S  CL lif-UTS  V  ’e 

eV  r*  r  0  r»C  r««,  or  r(3 


Figure  3.3  Energy  level  diagram  of  N  atom,  illustrated  the  same  as  Figure 
3.1.  (from  Vallance  Jones,  1974). 
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For  the  NI  1  series,  the  transitions  take  place  from  (l.v*’2.v 2/; )3/>  4D 
to  3.s  4P,  giving  a  total  of  8  allowed  transitions  or  emission  lines.  In  the  case 
of  NI  8,  the  transition  is  from  3/>  2P‘  to  3.v  P,  yielding  4  allowed  transitions 
or  emission  lines.  Table  3.1  gives  a  listing  of  each  wavelength  and  Einstein 
transition  probability  for  the  NI  1  and  NI  8  senes  of  emission  lines.  Notice 
only  the  strongest  lines  (the  ones  with  the  highest  Einstein  probability)  in 
the  NI  1  series  were  used  in  this  study  for  simplicity  sake.  It  was  also  felt 
that  the  weakest  lines  would  not  overcome  statistical  errors  to  be  analyzed 
in  a  consistent  manner.  Again,  ionic  emissions  were  not  analyzed. 

3.4.2  Analysis 

A  time  plot  of  each  of  the  nine  atomic  nitrogen  emission  lines  is 
shown  in  Figures  3.4  and  3.5.  It  is  immediately  obvious  that  all  of  the  lines 
follow  the  same  general  time  pattern.  The  only  exception  is  where  the  NI  1 
8680  and  8683  lines  change  possession  of  the  lead  near  18.9  EST  in  Figure 
3.4.  It  is  not  clear  if  this  is  an  anomaly  of  the  analysis  of  the  data  or  if  it  is 
indeed  truth. 

As  a  test  of  how  reasonable  the  analyzed  line  intensities  are,  we  can 
determine  the  ratio  of  two  lines,  and  compare  the  ratio  with  some 
experimental  results.  To  do  this  we  can  use  the  relation: 

i  =  3.6 

/,  77,  N,  A. 

where  r]  is  volume  emission  rate,  and  N  is  the  total  number  density  of  the 
upper  state. 
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LOCAL  TIME 


If  the  two  emissions  to  be  compared  are  from  the  same  upper  state, 
then  Ni  =  N2,  so  we  can  simplify  equation  3.6  to: 


which  is  just  the  ratio  of  the  Einstein  transition  probabilities,  and  we  do  not 
need  to  work  with  number  densities  at  all.  For  comparison  purposes,  two 
emission  lines  were  chosen  from  each  series,  NI  1  and  NI  8.  Using  the 
Einstein  probabilities  from  Table  3.1,  a  comparison  has  been  made  of  the 
ratios  from  the  two  sets  of  lines,  to  the  expected  values  obtained  from  the 
transition  probabilities.  Figures  3.6  and  3.7  show  the  results. 
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MEAN 


1  W 


.  •  /  Vl  .  .  5  .  f 

•  ,vr* f'  .  ..  ’  . ;  ,  „• 

■  -- _ _ _ _ _ _ _ — _ _ >.•  » _ .•  I- 1a.  ^  v.-.. 


•ACTUAL  RATIO 


18.5  1  9  19.5  20  20.5  21 

LOCAL  TIME 


Figure  3.6  Intensity  ratio  of  the  8594. 0A  to  8655. 9A  lines  of  the 
atomic  nitrogen  NI  8  series. 
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Figure  3.7  Intensity  ratio  of  the  8683.4A  to  8718.8A  lines  of  the  atomic 
nitrogen  NI  1  series. 


As  can  be  seen  from  Figures  3.6  and  3.7,  there  is  not  ideal  agreement 
between  the  analyzed  data  and  the  experimental  data,  taken  from  Wiese  et 
al.  The  analyzed  line  ratio  is  off  by  a  factor  of  1/2  in  Figure  3.6  and  a  factor 
of  about  3.5  in  Figure  3.7.  If  the  transition  probabilities  are  indeed 
accurate,  a  probable  explaination  for  these  deviations  is  that  the  synthetic 
spectrum  is  incompletely  analyzed,  resulting  in  improper  relative 
intensities.  This  may  be  partially  attributed  to  the  unexpectedly  large 
continuum  values  discussed  in  more  detail  in  chapters  4  and  5. 
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CHAPTER  4 


MOLECULAR  EMISSIONS 


4.1  Background 


4.1.1  Molecular  Term  Type  Notation 

Molecular  term  type  notation  is  very  similar  to  the  atomic  case 
discussed  in  chapter  3.  Both  cases  are  shown  here  for  comparison 
purposes: 


Atom 

2.V  +  1  , 


Lj 


Molecule 


2.V  +  I  / 


The  primary  difference  is  where  S,P,D,F  are  used  when  the  resultant 
angular  momentum  number,  L=0, 1,2,3  for  atoms,  molecular  notation  uses 
£,n,A,<J>  when  L=0, 1,2,3.  Just  as  for  atoms,  S  is  the  resultant  spin  quantum 
number,  and  Q.  (instead  of  J)  is  the  resultant  total  angular  momentum 
quantum  number  for  molecules.  Using  the  example  in  chapter  3  to  show 
the  similarity,  the  \P2I(,  would  be,  in  the  molecular  case,  3n2)(1. 
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4.1.2  Emission  Analysis 

An  advantage  of  analyzing  the  spectra  for  molecular  emissions  is 
that  one  can  recover  temperatures  from  the  rotational  distribution  of  each 
vibration-rotation  band  analyzed  without  a  high  resolution  instrument. 
This  requires  of  course  a  band  model  for  the  particular  type  of  transition 
involoved,  as  mentioned  in  chapter  2. 

The  molecular  emissions  analyzed  in  this  study,  in  order  of  intensity 
of  emission,  are:  N2  1st  Positive  (N2  IPG),  O2  Atmospheric  (O2  Atm),  N2 
Infrared  Afterglow  (N2  IRA),  N;  Meinel  ( N;  M),  and  hydroxyl  (OH).  The 
best  fit  obtained  from  the  31  parameters  listed  in  Table  2.3  is  shown  in 
figure  4.1. 

In  Figure  4.1,  the  highest  solid  line  is  the  corrected  measurement  at 
time  19.02  (19:01  EST),  the  dotted  line  near  the  measurement  line  is  the  total 
fitted  synthetic  spectrum  with  all  31  parameters'  contributions  added 
together,  though  only  the  molecular  contributions  are  shown  in  the  graph, 
located  below  the  data  and  fitted  lines.  The  atomic  species  were  discussed 
in  chapter  3. 

At  the  top  of  the  graph  are  the  band  origins  of  the  vibrational 
transitions  of  each  species  analyzed.  As  an  example  of  what  the  notation 
means,  the  bottom  line  shows  where  the  N2  IPG  (2,1)  band  origin  is  located. 

The  "2,1"  means  the  transition  takes  place  from  rotational  levels  of 
vibrational  level  2  of  the  upper  electronic  state,  to  rotational  levels  of 
vibrational  level  1  of  the  lower  electronic  state.  An  arrow  in  Figure  4.1 
indicates  band  origins  which  are  outside  the  range  of  the  data  but  still  have 
tails  which  contribute  to  the  overall  spectrum.  Figure  4.2  is  an  illustration 
of  all  three  types  of  energy  levels  in  a  molecule:  electronic,  vibrational,  and 
rotational. 
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Figure  4.1  Analyzed  molecular  composition  of  the  November  8,  1991  aurora 
over  Ann  Arbor,  MI.  The  top  solid  line  indicates  the  measurement, 
dotted  indicates  fitted  spectrum,  lower  lines  indicate  fitted  molecular 
parameters  and  continua. 


ionization,  excess  energy 
carried  away  by  photoelectron 


AB+  + 


AB 


inter-nuclear  distance 


Figure  4.2  Energy  level  of  a  typical  molecule  showing  electronic, 
vibrational,  and  rotational  energy  levels. 

The  solid  line  near  the  bottom  of  the  graph  is  the  continuum  used  to 
fit  the  spectrum's  "leftovers  ".  As  seen  in  Table  1.2  the  continuum  varied 
from  25.87  R/A  at  8400A  down  to  8.51  R/A  at  8900A.  The  25.87  R/A 
continuum  parameter  is  a  much  higher  value  than  expected.  Much  effort 
was  spent  trying  to  lower  this  parameter.  More  on  the  continuum  will  be 
discussed  in  chapter  5. 
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4.2  Nitrogen 


Diatomic  nitrogen  accounts  for  a  major  portion  of  auroral  emissions. 
Figure  4.3  shows  the  electronic  energy  states  of  N2.  Two  separate  N2,  and 
one  N2+,  band  systems  were  analyzed  in  the  8400-8900 A  spectrum:  1)  the 
first  positive  system,  which  occurs  due  to  B2Ug  -  AyI.*u  transitions  of  N2  as 

can  be  seen  in  Figure  4.3,  2)  the  infrared  afterglow  system  (IRA), 

occurring  due  to  B'3Z"  transitions,  and  3)  the  N2+  Meinel  system, 

from  A2  ri11(  -X2L*  transitions.  Each  system  will  be  covered  separately. 

4.2.1  First  Positive  System 

This  system  represents  the  most  intense  and  dominant  feature  in  the 
spectrum  considered.  The  production  of  the  B-state  is  produced  from 
energetic  electron  collisional  excitation  either  directly  from  the  X-state  (the 

ground  state)  or  indirectly  through  the  C  or  B'  states,  for  example, 
cascading  down  to  the  B-state  through  radiative  processes.  The 

transition  is  a  spin-allowed  electric  dipole  transition  and  occurs  very  fast: 
the  radiative  lifetime  of  the  B-state  is  on  the  order  of  105  seconds  (Loftus 
and  Krupenie,  1977),  and  the  primary  means  of  deexcitation  is  through  the 
first  positive  radiative  process. 
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eV 


Figure  4.3  Potential  energy  levels  and  electronic  band  systems  for  N2  and 

N2+  illustrated  the  same  as  Figure  3.1.  (From  Vallance  Jones,  1974). 

Five  parameters  were  used  in  this  study  to  account  for  the  system’s 
emissions:  four  band  brightnesses  (the  (1,0), (2,1), (3, 2),  and  (4,3)  bands)  and 
one  temperature  which  would  be  used  for  all  four  bands.  As  can  be  seen 
from  Table  2.3,  the  fit  of  the  (4,3)  band  was  far  from  ideal  (noted  by  the 
error).  This  could  be  due  to  a  number  of  reasons:  1)  low  band  brightness,  2) 
band  origin  near  the  limit  of  the  measured  spectrum,  or  3)  the  abnormally 
high  continuum  at  that  end  "stealling"  some  of  the  intensity  from  the  (4,3) 
band.  The  band  was  left  in  the  analysis  because  it  was  felt  it  would  help 
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keep  the  continuum  lower  in  magnitude.  In  any  event,  the  poor  fit  of  the 
(4,3)  band  requires  leaving  it  out  of  the  graphs  that  follow  for  clarity's  sake. 

Figure  4.4  shows  a  summary  of  the  analyzed  parameters.  The  scale 
for  the  three  band  brightnesses  is  on  the  left,  while  the  scale  for  the 
temperature  is  on  the  right.  Though  it  must  be  stressed  that  some  of  the 
variations  in  brightnesses  and  temperatures  may  be  due  to  analytical 
fluctuations  versus  physical,  general  observations  may  be  taken  from  the 
graphs.  Quite  possibly  most  of  the  variations  may  be  later  confirmed, 
through  comparisons  with  other  data,  as  real. 

The  most  obvious  feature  of  Figure  4.4  is  that  all  three  brightnesses 
follow  nearly  the  exact  same  pattern.  This  is  comforting  and  suggests  a 
somewhat  constant  vibrational  distribution  with  time.  This  is  further 
investigated  later  in  this  section. 

As  seen  with  the  atomic  emissions,  the  pattern  of  brightnesses  in 
Figure  4.4  shows  about  a  25  minute  periodicity.  It  can  also  be  seen  that  the 
temperature  variation  does  not  correspond  to  brightness  variation.  In  fact 
temperature  and  brightness  seem  to  be  anticorrelated,  noting  their 
relationships  at  times  19.0  and  19.9  in  the  figure,  in  particular. 

As  mentioned  shortly  before,  the  band  brightnesses  varied  similarly 
with  respect  to  one  another  over  time.  This  fairly  constant  variation  allows 
us  to  make  a  vibrational  population  comparison.  Cartwright  et  al.  (1971) 
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Figure  4.4  N2  first  positive  band  brightnesses  and  temperature 
recovered  from  the  aurora  measurements. 


made  a  comparison  of  the  vibrational  population  model  they  theorized  with 
that  of  Shemansky.  Here  I  have  simply  extracted  the  data  from  both  models 
up  to  v=6,  and  superimposed  the  data  derived  from  this  study  averaged  over 
the  period  from  18.5  to  about  21.0  EST  for  comparison  purposes.  Figure  4.5 
shows  the  results. 


Figure  4.5  N2  A-state  vibrational  population  derived  from  first  positive 
band  brightness  ratios.  Populations  are  normalized  to  the  v=3 
level.  Shemansky  and  Cartwright  data  taken  from  Cartwright 
et  al.,  1971. 


As  can  be  seen  from  the  graph,  the  v=l,2  populations  from  this  study 
seem  to  be  too  high  when  normalized  to  the  v=3  level.  This  means  that 
either  the  (1,0)  and  (2,1)  bands  are  too  bright,  according  to  the  two  theories, 
or  the  (3,2)  band  is  too  dim.  I  believe  the  latter. 

When  viewing  Figure  4.1,  one  can  see  that  the  high  continuum 
mentioned  earlier  slopes  from  high  to  low  wavelengths.  Since  the  (2,1)  and 
(1,0)  bands  lie  to  the  high  wavelength  side  of  the  (3,2)  band,  their  band 
brightnesses  are  not  affected  as  much  as  the  (3,2)  band  by  the  sloped 
continuum.  I  believe  that  were  it  not  for  the  unexpectedly  high  amplitude 
of  the  continuum,  the  ratios  of  the  intensities  of  the  (1,0)  and  (2,1)  to  (3,2) 


bands  would  be  much  closer  to  the  values  predicted  by  Cartwright,  in  a 
qualitative  sense,  at  least.  However  with  the  uncertainty  of  fit  of  the  (3,2) 
band  at  4%,  and  considering  the  slope  of  the  data  in  Figure  4.5,  agreement 
with  Shemansky's  vibrational  distribution  is  not  out  of  the  question. 

4.2.2  Infrared  Afterglow 

This  system  is  one  of  the  newest  of  the  emissions  known  to  occur  in 
the  aurora.  Due  primarily  to  overshadowing  by  first  positive  emissions,  the 
infrared  afterglow  system,  and  consequently  the  upper  state  (B  :JXU  ■),  were 

not  known  about  until  the  1950's.  The  molecular  constants  required  for 
band  models  to  predict  these  emissions,  were  not  available  until  around 
1960.  Since  then,  the  N2  IRA  emissions  have  been  analyzed  numerous 

times  (Gattinger  and  Vallance  Jones,  1974,  for  example). 

A  view  of  the  states  involved  in  the  transition  ( B' '  X(,  -B‘ I  I,  )  can  be 

seen  in  Figure  4.3.  Excitation  of  the  N2  molecule  to  the  B -state  is  primarily 
through  energetic  electron  excitation.  There  are  few  radiative  transitions 
to  the  B'-state,  so  cascading  to  this  state  should  not  occur  often. 

A  large  portion  of  this  study  was  spent  trying  to  decrease  the 
magnitude  of  the  continuum  and  to  lower  the  value  of  y2.  I  believed  the  N2 
IRA  system  would  resolve  most  of  the  continuum  magnitude  and  high  y2 
problems.  Much  time  was  spent  trying  to  develop  a  band  model  for  the 
system.  None  could  be  located  in  time  for  use,  when  it  was  found  that  it 
would  be  possible  to  build  such  a  model  quite  easily,  noting  that  the  <1  -  M  I 
transition  is  just  the  onposite  of  the  first  positive's:  TI  -  'X. 

Since  E.  Buscela  (a  PhD  candidate  at  the  University  of  Michigan) 
already  had  a  model  for  the  first  positive  system  which  was  easily  alterable, 
his  model  was  chosen  to  be  used  for  the  origin  of  an  N2  IRA  model.  All  that 


was  done  to  make  the  switch  from  first  positive  to  IRA  was  to  change  the 
designation  of  the  line  strengths  so  that,  for  example,  a  rotational 
transition  from  level  2  to  1  of  the  '□  -  :,Z  became  a  1  to  2  rotational  level 
transition  of  the  :tI  -  ‘fl.  Therefore  a  P-branch  of  the  first  positive  system 
became  an  R-branch  of  the  IRA  system,  Q-branch  line  strength  indeces 
stayed  the  same,  and  l-branches  became  P-branches.  A  test  of  the  new 
model  was  accomplished  and  Figure  4.6  shows  the  comparison. 

Figure  4.6  was  created  by  convolving  the  individual  rotational  lines 
into  the  NIRS  instrument  function  so  that  a  channel  *  channel 
comparison  of  model  vs  data  could  be  made.  Then  a  non-linear  least 
squares  fit  of  the  band  model  to  the  data  for  both  intensity  and  temperature 
was  done.  A  band  brightness  of  1280  R  and  temperature  of  208  K  were 
found  to  give  the  model  the  best  fit  to  the  data.  Molecular  constants  were 
obtained  from  Huber  and  Herzberg  ( 1979). 

It  is  evident  that  the  model  does  a  decent  job  of  predicting  the  spectral 
shape  and  positions  of  the  IRA  bands,  though  I  am  not  fully  satisfied.  The 
P,  Q,  and  R-branch  peaks  appear  to  be  displaced  in  wavelength  somewhat 
(about  4  A),  and  the  intensity  of  the  peaks  are  not  in  proper  proportion  with 
one  another.  In  addition,  the  results  of  the  analysis  using  the  three  IRA 
bands  shown  in  Table  2.3  indicates  the  analytical  error  of  the  IRA  bands  is 
about  twice  that  of  other  bands  of  similar  brightnesses,  and  the  consistency 
of  the  IRA  recovered  parameters  were  mucn  less  consistent  than  other 
species’  parameters. 

Nonetheless,  the  IRA  model  did  improve  the  fit  of  the  synthetic 
spectrum  to  the  measurements  somewhat,  so  the  model  was  used  as  is,  but 
with  caution.  The  IRA  band  brightnesses  and  temperature  recovered  from 
the  analysis  were  found  to  be  too  eratic  to  be  shown  in  graphic  format. 
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CONVOLVED  N,  IRA  (8,3)  BAND 


>ire  4.6  Comparison  of  modelled  N2  IRA  with  experimental  data  of 
Dieke  and  Heath  (1960). 


4.2.3  N.,+  Meinel 


This  system  was  the  only  ion  species  analyzed  in  this  study.  It  has 
been  widely  studied  in  the  aurora  and  known  to  exist  for  quite  some  time. 
The  primary  source  of  production  of  the  upper  state  of  the  Meinel  band 
transition  (A2nii(-X2E*)  is  electron  impact  ionization  excitation  (Rees, 

1987): 


e  +  N,  — >  Nj(A  2nu)  +  e.  4.1 

For  this  system,  only  two  parameters  were  chosen  to  be  analyzed:  the 
(5,3)  and  (4,2)  band  brightnesses.  Since  our  spectral  measurements  do  not 
have  the  sensitivity  and  resolution  to  produce  a  meaningful  temperature, 
the  temperature  here  was  fixed  at  600K. 

Figure  4.7  shows  the  variation  of  the  (4,2)  and  (5,3)  band  brightnesses 
with  time.  Although  this  graph  has  a  linear  scale  on  the  vertical  and  the 
previous  graphs  have  a  log  scale,  one  can  still  compare  this  graph  with  the 
others  in  a  general  sense  and  see  the  brightnesses  follow  the  same  basic 
pattern. 

4.3  Oxygen 

The  only  O2  band  system  analyzed  in  this  study  is  that  of  the 
atmospheric  system.  Although  the  /?'!*- X'Z~  transition  is  a  spin 

forbidden,  magnetic  dipole  process,  its  appearance  in  the  aurora  is  quite 
prominant.  Observations  of  bands  of  this  system  were  first  taken  from 
solar  measurents  back  in  the  early  1800’s  (Krupenie,  1972).  The  energy 
levels  involved  in  the  transition  are  shown  in  Figure  4.8. 


18.5  19.0  19.5  20.0  20.5  21.0  21.5 

LOCAL  TIME 


Figure  4.7  Time  series  of  N2+  Meinel  band  brightnesses  recovered  from 
the  aurora  measurements. 


The  mechanism  which  produces  the  excited  b-state  of  O2  in  the  lower 
thermosphere  is  a  chemical  process.  During  nighttime  and  magnetically 
quiet  periods,  the  emission  occurs  near  95  km  and  is  produced  from: 

0  +  0  + m  ->o;  +  M 

4  2 

o;+o, -»o2(b ‘z;)+o2. 

At  other  times,  such  as  in  the  dayglow  and  and  the  aurora,  the  production 
of  02(b1L)  is  preceded  by  the  production  of  O^D): 


GO 


O  +  e— »  0(‘D)  +  e 

O; +e -><)(’ D)  +  ()  4.3 

O,  +/jv->0(‘D)  +  0, 

which  leads  to: 

0(,D)  +  03  ->02(h,I)  +  0(,P),  4.4 

which  then  can  radiate  in  the  atmospheric  system.  Equation  4.4  is  the 
main  source  of  these  emissions  in  the  aurora. 


Figure  4.8  Potential  energy  levels  and  electronic  band  systems  for  O2  and 
02+  illustrated  the  same  as  Figure  3.1.  (From  Vallance  Jones,  1974). 
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The  lifetime  of  the  state  is  on  the  order  of  about  10  seconds.  Since 
this  state  has  such  a  long  lifetime,  many  collisions  with  the  ambient  gas 
will  take  place  before  the  emission  occurs.  We  can  then  say  the  O2 
atmospheric  temperatures  recovered  from  the  measurement  will  be  pretty 
representative  of  the  kinetic  temperature  of  the  ambient  gas  where  the 
emission  is  taking  place.  This  is  not  necessarily  true  for  the  recovered  N2 
temperatures,  which  primarily  represent  the  rotational  distribution,  owing 
to  the  short  lifetimes  of  the  upper  states,  but  we  can  still  say  the  N2 
temperatures  approximate  the  kinetic  temperature. 

Figure  4.9  illustrates  the  time  variation  of  the  recovered  O2 
atmospheric  parameters.  Three  band  brightnesses  and  two  temperatures 
were  chosen.  Two  temperatures  were  used  in  this  case  to  see  if  more  than 
one  temperature  could  be  recovered  successfully,  and  indeed  it  could  be 
done.  This  was  not  attempted  with  all  other  bands  to  keep  computation 
time  minimal. 

The  overall  shape  of  the  brightnesses  and  temperatures  compare 
well  among  themselves  and  with  those  of  the  N2  parameters  as  well.  It  is 
evident  that  the  larger  errors  present  in  the  weaker  (2,1)  and  (3,2)  bands 
creates  slightly  greater  variability  in  both  temperature  and  brightnesses. 

It  is  interesting  to  note  the  relationship  of  the  two  temperature  plots 
of  Figure  4.9.  The  average  temperature  of  the  (0,1)  band  is  400K,  while  the 
average  of  the  (1,2)  and  (2,3)  bands’  combined  temperature  is  625K.  This  is 
as  we  might  expect.  The  higher  vibrational  states  are  more  likely  to  exist  at 
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Figure  4.9  Time  series  of  02  atmospheric  band  brightnesses  and 
temperatures  recovered  from  the  aurora  measurements. 

higher  altitudes  in  the  thermosphere  where  there  are  fewer  collisions  that 
would  deexcite,  or  quench,  those  states.  We  also  know  the  ambient 
temperature  increases  with  height  in  the  thermosphere,  where  auroral 
emissions  like  the  O2  atmospheric  emissions  originate.  We  would  therefore 
expect  the  temperature  of  the  (1,2)  and  (2,3)  bands  to  be  higher  than  the  (0,1) 
band  temperature. 
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4.4  Hydroxyl 


The  hydroxyl  radical's  (chemical  abbreviation  OH)  emissions  are  a 
very  widespread  and  dominant  feature  in  any  nightglow  measurement  of 
the  infrared  spectrum.  Unlike  the  homonuclear  molecules,  OH  has  a 
permanent  dipole  moment  that  allows  the  molecule  to  make  vibration- 
rotation  as  well  as  electronic-vibration-rotation  transitions.  Thus  it  can 
have  more  band  emissions  than  diatomic  molecules  in  the  infrared 
spectrum. 

The  emissions,  in  the  portion  of  the  spectrum  concerned  with  here, 

are  due  to  overtone  (Ao>l)  transitions  from  vibrationally  excited  radicals  in 
the  ground  electronic  state  of  OH:  X2n,,.  The  primary  source  of  the 

vibrationally  excited  radicals  is  the  displacement  reaction  between  H  and  O3 
(Bates  and  Nicolet,  1950): 

H  +  O,  -»  0H(v’<9)  +  0,  +3.3eV  .  4.4 

The  3.3  eV  of  energy  released  by  the  reaction  is  enough  to  excite  the  radical 
up  to  and  including  the  v'= 9  vibrational  level,  sufficient  to  excite  the 
molecule  to  the  v’=6,7  levels  investigated  here. 

Altitude  profiles  of  the  Meinel  bands  obtained  by  rocket  borne 
instruments  place  the  peak  emission  between  about  85  and  90  km  altitude. 
Since  the  rotational  relaxation  of  OH  is  believed  to  be  rapid  enough  to 
produce  a  thermal  distribution  of  rotational  populations,  analysis  of  the 
rotational  distribution  yields  the  kinetic  temperature  of  the  ambient  gas,  as 
with  the  O2  atmospheric  system  (Rees,  1989).  Because  the  peak  emission 
altitude  of  the  OH  Meinel  bands  is  below  the  thermosphere,  and  the 


emissions  indicate  ambient  conditions  of  the  5  km  altitude  layer  from 
where  the  photons  originate,  the  OH  Meinel  bands  are  used  extensively  for 
investigating  the  effects  of  gravity  wave  propagation  through  the  region. 

Figure  4.10  shows  the  OH  parameters  recovered  from  the  aurora 
measurements.  It  is  immediately  obvious  that  these  emissons  do  not  follow 
the  same  pattern  of  any  of  the  species  shown  up  to  this  point.  This  is  as 
expected:  since  the  peak  emission  altitude  of  the  OH  Meinel  bands  is  near 
85  km,  and  auroral  electrons  lose  the  bulk  of  their  energy  by  inelastic 
scattering  above  100  km,  we  would  not  expect  OH  Meinel  emissions  to 
immediately  follow  the  others,  which  emanate  from  above  the  mesopause. 
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Figure  4.10  Time  series  of  OH  band  temperature  and  brightnesses 
recovered  from  the  aurora  measurements. 


4.6  Comparisons 


4.6.1  Temperatures 

One  particularly  useful  indication  temperature  gives  us  is  the 
altitude  from  which  the  emission  is  taking  place.  By  using  a  model  such  as 
the  empirical  MSIS-86  model  (Hedin,  1987),  we  can  predict  the  temperature 
profile  of  the  lower  thermosphere.  This  will  then  allow  us  to  make  an 
estimation  on  the  emission  altitude  by  matching  the  temperature  with  a 
temperature  profile  from  MSIS-86. 

Figure  4.11  shows  a  plot  of  the  temperature  profile  output  from  MSIS- 
86.  The  parameters  input  to  the  model  to  get  this  plot  were:  date=91313, 
time=0.0  UT  (8  PM  EST),  F10.7=200,  F10.7A=197,  Ap=65,  latitude=42, 
longitude=85  (Ann  Arbor).  Thus  the  profile  is  a  snapshot,  over  Ann  Arbor, 
MI,  of  the  lower  thermosphere  at  8  PM  on  the  evening  of  Nov.  8,  the  center 
in  time  of  the  aurora  measurements. 

Also  included  in  Figure  4.11  are  the  four  average  temperatures 
during  the  aurora,  recovered  from  the  analysis,  shown  in  previous  plots. 
The  N2  IRA  temperature  again  has  been  left  out  because  of  its  high 
variability.  The  four  temperatures  have  been  assigned  to  an  altitude  on  the 
MSIS-86  profile  so  that  we  can  estimate  where  the  peak  of  the  emission 
layer  is  occurring. 


TEMPERATURE  PROFILE 


Figure  4.11  Temperature  profile  from  85  -  200  km.  Profile  is  from  the 
MSIS-86  empirical  model,  highlighted  points  on  the  graph  are 
average  species  temperatures  during  the  aurora. 

Figure  4.12  is  related  to  Figure  4.11,  showing  the  variation  of  the 
recovered  temperatures  with  time,  with  specific  altitude  (95,115,125,135,150) 
temperatures,  as  a  function  of  time,  added  to  show  how  the  peak  emission 
altitude  of  each  species  varied  in  that  time  period.  The  five  altitudes  levels 
were  chosen  simply  on  the  basis  of  spacing. 

If  we  assume  that  the  changes  in  temperature  of  the  species  are 
primarily  due  to  changes  in  the  altitude  of  emission,  we  can  see  how  the 
peak  emission  altitudes  changed  with  time  in  Figure  4.12.  The  variation 
for  the  three  thermospheric  emissions  is  on  the  order  of  20-30  kilometers, 
while  the  variation  of  OH's  peak  emission  altitude  is  less  than  10  kilometers 
between  maximum  and  minimum  peaks. 
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Figure  4.12  Comparison  of  species  temperatures  and  MSIS-86  altitude 
temperatures. 


Since  the  hydroxyl  radical's  emission  is  largely  unaffected  by  the 
auroral  electrons,  it's  temperature  variation  is  probably  primarily  due  to 
actual  changes  in  kinetic  temperature  and  secondarily  to  emission  altitude. 
The  OH  temperature  variation  is  considered  to  be  dominated  by  upward 
propagating  gravity  waves,  and  a  spectral  analysis  on  the  temperature 
variation  would  yield  useful  information  on  the  influence  of  the  waves. 

One  last  note  on  the  hydroxyl  emission:  as  was  noted  earlier,  the 
peak  emission  altitude  of  OH  is  generally  considered  to  be  between  85  and  90 


km,  as  measured  by  rocket-borne  instruments.  The  temperature  recovered 
from  the  aurora  measurements  averaged  to  about  180K,  but  the  MSIS-86 
model  predicted  a  temperature  at  90  km  to  be  196  K,  and  a  minimum 
temperature  at  the  mesopause  (95  km)  to  be  186  K.  Thus  the  recovered 
temperatures  from  the  OH  bands  is  lower  than  any  altitude  s  temperature 
predicted  by  MSIS.  It  is  known  that  the  model’s  accuracy  is  +  or  -  10-15% 
(about  3  K).  The  recovered  temperature  error  is  probably  good  to  within 
about  5  K.  With  these  error  limits  in  mind,  the  recovered  temperatures  still 
suggest  a  peak  emission  altitude  greater  than  90  km. 

4.6.2  Temperature  vs.  Brightness 

Interesting  observations  can  be  made  when  comparing  species 
brightness  with  it's  temperature,  as  done  before,  but  with  a  little  more 
detail.  Figures  4.13  ana  4.14  show  band  brightnesses  versus  species 
temperature  for  O2  and  N2  respectively. 

We  noted  before  that  the  temperatures  and  brightnesses  seemed  to  be 
anticorrelated,  or  out  of  phase  by  about  180  degrees.  From  the  figure  above, 
with  the  scales  expanded,  one  can  see  this  is  definitely  the  case.  The  two 
plots  are  not  perfectly  out  of  phase,  but  close  to  it  over  most  of  the  period, 
within  about  three  minutes. 

Our  results  suggest  that  when  emissions  are  their  brightest,  which 
is  when  the  auroral  electron  flux  is  greatest,  the  ambient  temperature  is 
lowest.  This  notion  can  be  explained  in  the  following  way.  A  burst  of 
highly  energetic  auroral  electrons  precipitating  downward  can  penetrate 
deeper  into  the  atmosphere  than  a  burst  of  less  energetic  electrons.  This  is 
because  more  highly  energetic  electrons  can  experience  more  inelastic 
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Figure  4.13  O2  Atmospheric  (0,1)  band  brightness  versus  temperature. 
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Figure  4.14  N,>  first  positive  ( 1 ,0)  band  brightness  versus  temperature 
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collisions  with  ambient  particles  on  their  way  downward  through  the 
earth's  atmosphere,  before  giving  up  all  their  excess  kinetic  energy,  than 
less  energetic  particles  can.  So  the  greater  the  average  energy  of  an 
auroral  electron  distribution,  the  deeper  the  energy  can  penetrate,  and  the 
peak  emission  will  come  from  a  lower  altitude  due  to  the  greater  density  at 
lower  altitudes.  In  addition,  the  brightness  of  the  aurora  is  determined  by 
the  magnitude  of  the  flux  of  auroral  electrons  (expressed  in  terms  of 
ergs/cm2/sec  in  cgs  units).  Since  the  recovered  data  suggests  lower 
temperatures  (and  therefore  emission  altitude)  accompany  higher 
intensities  in  the  aurora,  then  the  fluxes  of  more  energetic  electrons  must 
have  come  in  greater  numbers  also. 

To  summarize,  1)  the  temperature  suggests  a  penetration  depth  of 
electrons.  2)  The  penetration  depth  suggests  the  characteristic  energy  of 
the  auroral  particles.  3)  The  magnitude  of  the  flux  of  auroral  electrons 
affects  the  brightness  of  the  aurora.  One  can  therefore  use  the  recovered 
temperatures  and  brightnesses  to  estimate  the  characteristic  energy  and 
flux  of  auroral  particles  through  the  use  of  a  model. 

This  concludes  chapter  4  and  the  entire  study's  analysis  of  the 
aurora  measurements.  The  final  chapter,  chapter  5,  will  restate  what  I 
have  already  presented,  and  will  suggest  areas  where  improvement  can  be 
made  so  that  an  even  more  complete  analysis  of  the  data,  and  verification  of 
resuUs  can  be  achieved. 
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CHAPTER  5 


CONCLUSIONS 


5.1  General  Remarks 

On  the  evening  of  November  8,  1991,  a  moderate  intensity  auroral 
event  was  observed  over  Ann  Arbor,  MI  (42  N,  85  W),  as  well  as  many  more 
equatorward  locations.  Since  this  region  is  outside  the  auroral  zone,  the 
event  is  a  fairly  uncommon  one  at  this  location,  and  if  the  measurements 
had  been  taken  inside  the  auroral  zone,  the  intensities  would  have  been 
extremely  bright,  much  more  so  than  measured  in  Ann  Arbor.  This  fairly 
unusual  event  for  this  region  was  measured  by  the  University  of 
Michigan's  Near-Infrared  Spectrometer  from  6:30  to  12:00  pm  EST 
continuously  for  integration  periods  of  30  seconds.  A  description  of  the 
instrument  and  it's  CCD  detector  was  given  in  chapter  1 . 

The  instrument  resolution  is  about  13  A,  and  NIRS  was  set  up  to 
measure  the  spectral  region  of  approximately  8400-8900  A  for  a  nightglow 
study  when  the  unexpected  auroral  event  occurred.  With  a  30  second 
integration  time,  and  considering  only  photon  statistics,  the  measurement 
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readings  at  each  of  the  384  data  channels  are  accurate  to  within  +/-  0.5%  for 
brighter  periods  and  +/-  1.5%  for  less  intense  periods. 

Retrieval  of  useful  emission  parameters  such  as  atomic  line 
intensities,  molecular  band  brightness,  and  molecular  band  temperature 
require  analysis  of  the  data  and  the  creation  of  a  synthetic  spectrum,  which 
accounts  for  all  emissions  believed  to  exist  in  the  measurement.  The 
inclusion  of  molecular  emissions  in  the  synthetic  spectrum  requires  a  band 
model  for  each  type  of  transition.  The  band  model  output  hundreds  of 
individual  rotational  lines  emitted  due  to  the  molecule’s  transition,  based 
on  the  inputted  total  band  brightness  and  temperature. 

Once  the  entire  synthetic  spectrum  is  constructed  from  atomic  and 
molecular  parameters,  then  convolved,  analysis  of  the  data  through  the  use 
of  a  least  squares  technique  is  required,  so  that  the  best  possible  fit  of  the 
synthetic  spectrum  to  the  measurement  is  found.  A  non-linear  least 
squares  technique  discussed  in  chapter  two  was  used  for  this  study  to 
recover  a  total  of  31  parameters  from  the  measurements,  including  low-end 
and  high-end  continuum  brightness  parameters  used  to  "fill  in  the  gaps". 

The  iterative  non-linear  least  squares  technique  was  used  on  the 
entire  evening's  measurements,  about  600  30-second  frames.  The  vast 
majority  of  this  study  concentrated  on  the  time  period  from  6:30  to  9:30  pm 
EST,  when  the  aurora  was  at  it's  brightest.  Once  I  was  satisfied  with  the 
goodness  of  fit  of  the  synthetic  spectrum  to  the  measurements,  judged 
primarily  from  the  value  of  chi-squared,  investigation  of  the  recovered 
parameters  for  physical  insight  was  begun. 

The  primary  purposes  of  this  study,  then,  was  to  build  a  non-linear 
least  squares  fitting  technique  that  was  flexible  enough  to  change  allow 
changes  in  parameters  with  minimal  effort,  arrive  at  a  suitable  synthetic 
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spectrum,  and  attempt  to  verify  physically  what  the  recovered  parameters 
implied. 

5.2  Results 

Chapter  3  discussed  features  found  in  the  atomic  species’  emissions. 
One  oxygen  and  nine  nitrogen  lines,  for  a  total  of  ten,  atomic  emission  lines 
were  recovei'ed  from  the  measurements.  The  most  dominant  lines  in  the 

o 

8400-8900A  spectral  range  measured,  in  order  of  most  to  least  intense, 
were:  01  4  8446A,  NI  1  8683,  NI  1  8680.  These  and  the  seven  other  atomic 
lines  were  plotted  as  a  function  of  time.  All  the  lines  showed  a  similar 
intensity  versus  time  pattern,  with  a  periodicity  of  about  25  minutes,  which 
was  apparently  driven  by  the  variation  of  the  flux  of  auroral  electrons. 

An  attempt  was  made  to  verify  the  relative  strengths  of  some  of  the 
nitrogen  lines  measured  To  do  so,  the  ratio  of  lines  with  the  same 
multiplicity  within  the  same  electronic  transition,  were  compared  with  the 
ratio  of  the  Einstein  transition  probabilities  (from  Wiese,  et  al.,  1966).  The 
results  showed  the  8594/8656A  intensity  ratio  to  be  lower  than  the  expected 
value  by  a  factor  of  about  1/2,  and  the  8683/87 19A  intensity  ratio  to  be  too 
high  by  a  factor  of  3.5. 

Chapter  4  discussed  the  molecular  species'  parameters  recovered 
from  the  aurora  data.  A  total  of  14  molecular  band  brightnesses  and  5  band 
temperatures  were  recovered  from  the  data.  The  molecular  species 
analyzed  in  this  study,  in  order  of  total  intensity  of  emission,  are:  Nz  1st 
Positive  (N2  IPG),  O2  Atmospheric  (O2  Atm),  N2  Infrared  Afterglow  (N2 
IRA),  N)  Meinel  ( M),  and  the  hydroxyl  radical  (OH  ). 


The  N2  IPG  system  was  the  most  useful  system  of  all  to  this  study, 
because  of  the  number  of  vibration-rotation  bands  (the  (1,0),  (2,1),  (3,2),  and 
(4,3)  bands)  in  the  measured  spectrum.  The  (4,3)  band  was  evidently  too 
weak,  or  too  far  off  of  the  low  wavelength  side  of  the  spectrum  to  be 
recovered  with  any  reliability.  The  other  three  bands,  however,  were 
consistent  enough  to  allow  a  study  of  vibrational  population  to  be 
accomplished,  using  data  taken  from  Cartwright,  et  al.  (1971). 

Our  results  gave  higher  populations  for  the  v=l,2  levels,  when 
normalized  to  v=3,  than  either  Cartwright's  or  Shemansky’s  modeled 
vibrational  population  in  the  aurora,  by  about  10%  and  20%,  respectively. 
While  the  values  of  the  populations  agreed  with  Cartwright's  model  better 
than  Shemansky's,  the  slope  of  this  study’s  vibrational  population, 
however,  tended  to  agree  with  the  model  of  Shewmansky  better.  It  was 
suggested  that  some  of  the  disagreement  of  the  results  with  both  models 
may  have  been  due  to  the  unexpectedly  high  value  of  the  recovered 
continuum  at  the  low  wavelength  side  of  the  spectrum. 

The  N2  IRA  system  was  added  to  the  list  of  parameters  when  it  was 

noted  how  many  bands  this  system  had  in  the  spectrum.  A  total  of  six 
bands  exist  in  or  near  the  8400-8900A  range.  It  was  thought  that  some,  if 
not  all,  of  these  bands  would  lower  the  high  continuum  values  and  also 
enhance  the  goodness  of  fit  by  lowering  x2-  After  trying  to  fit  all  the  bands, 
with  many  different  combinations,  it  was  found  that  the  system  only 
slightly  helped  lower  the  continuum  values  and  x2  In  addition,  the  error  of 
fitting  and  consistency  of  the  recovered  IRA  parameters  prohibited 
meaningful  examination  of  the  parameters.  It  was  felt  that  improvement 
of  the  band  model  developed  may  improve  conditions. 
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The  N2+  Meinel  system  has  two  bands  which  emit  in  the  spectrum 
studied.  These  two  bands  were  analyzed  with  a  fixed  temperature  of  600K. 
A  time  series  of  the  band  intensities  was  presented.  It  showed  a  similar 
band  brightness  versus  time  pattern  as  the  other  species. 

The  02  atmospheric  system's  three  bands  were  fitted  using  two 

temperatures:  one  corresponding  to  the  (0,1)  band  only,  and  one 
corresponding  to  the  (1,2)  and  (2,3)  bands.  The  results  give  a  higher 
temperature  for  the  (l,2)/(2,3)  bands  than  the  (0,1)  band.  This  suggests  a 
higher  peak  emission  altitude  for  the  higher  vibrational  states  (v’=l,2)  than 
the  v'=0  vibrational  state. 

Two  bands  of  the  hydroxyl  radical's  (OH)  Meinel  system  were 
analyzed  from  the  measi  "ements.  The  peak  emission  altitude  of  OH  (85-90 
km)  precludes  these  emissions  from  being  enhanced  by  auroral  electrons. 
The  recovered  band  brightnesses  and  temperatures  did  not  vary  in  a  way 
which  was  similar  to  that  of  species  associated  with  the  aurora. 

Comparisons  of  different  species'  temperatures  were  made  with 
reference  to  the  MSIS-86  temperature  profile.  The  average  altitude  of  peak 
emission  for  each  species  was:  OH-95  km,  N2  IPG-  115  km,  02  atm  (0,1)- 
117  km,  and  02  atm  (l,2)/(2,3)-  135  km.  The  OH  emission  altitude,  in 
particular,  is  notably  higher  than  anticipated. 

A  comparison  of  species  temperatures  with  MSIS-86  altitude 
temperatures  was  made,  in  order  to  show  variation  of  peak  emission 
altitude,  assuming  there  is  little  change  in  kinetic  temperature  taking 
place.  Results  suggest  variations  of  emission  altitude  are  about  25  km  for 
most  species,  while  the  variaton  of  hydroxyl  is  less  than  10  kilometers 
between  maximum  and  minimum  values. 
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When  closely  comparing  species  brightness  and  temperature  versus 
time,  one  can  see  the  recovered  data  shows  that  temperature  is  out  of  phase 
with  band  brightness.  Since,  from  the  discussion  in  section  4.6.2,  we  know 
that  species  temperature  can  yield  the  characteristic  energy  of  auroral 
particles,  and  the  magnitude  of  the  flux  of  auroral  electrons  determines 
brightness.  The  out  of  phase  relationship  between  temperature  and 
brightness  suggests  a  higher  characteristic  energy  of  auroral  electrons 
accompanied  greater  periods  of  flux  of  electrons  during  these 
measurements. 

5.2  Areas  for  Improvement 

As  suggested  in  various  places  throughout  this  writing,  there  are  a 
number  of  places  where  improvement  can  be  made  to  more  fully  perfect 
recovery  of  emission  species  parameters,  and  to  carry  further  the  use  of 
recovered  information  to  further  ends.  The  most  immediate  of  which  is  the 
modelling  of  the  auroral  particle  flux  using  temperatures  and  brightnesses 
recovered  from  the  data.  This,  when  combined  with  other  measurements, 
can  lead  to  a  greater  understanding  of  cause  and  effect  of  the  many 
variables  associated  with  the  aurora. 

Improvements  on  the  fitting  of  parameters  can  be  made.  As  noted 
earlier,  much  time  was  devoted  to  attempting  to  lower  y2  and  continuum 
magnitudes.  I  believe  that  most  of  the  problems  with  the  fitting  is  on  the 
low  wavelength  side,  due  to  the  magnitude  of  continuum  being  so  high  on 
that  end,  and  because  of  the  poor  fit  of  the  synthetic  spectrum  to  the 
measurement  near  8430A,  as  depicted  in  Figure  4.1.  There  are  two  possible 
explanations  for  the  less  than  idea!  fit:  1)  there  is  one  or  more  atomic  and/or 


78 


molecular  species  which  are  important  here  that  have  not  been  included  in 
the  analysis,  2)  there  may  be  some  internal  scattering  of  radiation  in  the 
NIRS  spectrometer  which  may  be  creating  "white  noise"  in  the 
measurements. 

The  first  option,  that  there  may  be  species  left  out  of  the  analysis,  is  a 
tough  one  to  solve.  Possible  species  which  have  not  been  tried,  due  to  the 
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belief  that  they  would  have  a  very  minor  effect,  are:  the  01  4  8426  A  atomic 
line,  and  the  (0,5)  band  of  the  02+  1NG  (first  negative  ground)  system.  The 

other  11-  NI  1  and  Nil  atomic  lines  may  improve  fitting  at  other  parts  of  the 
spectrum. 

The  possibility  of  white  noise,  or  internal  scattering  of  light  in  the 
spectrometer,  can  be  partly  justified  when  examining  Figure  5.1.  It  shows 
how  x2,  a  measure  of  goodness  of  fit,  varies  with  time.  An  ideal  fit  would 
have  a  x2  value  of  1.  The  01  8446  line  is  added  as  a  reference,  to  show  the 
intensity  of  the  aurora  with  time.  One  can  see  that  when  the  aurora 
brightens,  as  indicated  by  the  8446  line,  the  value  of  %2  goes  up  as  well,  and 
during  the  dimmer  periods  the  value  of  x2  goes  down. 

It  is  obvious  that  during  brighter  periods  there  will  be  more 
scattering  of  white  light  in  the  spectrometer,  and  during  less  bright  periods 
there  will  be  less  scattering,  if  it  is  indeed  occurring.  Since  the  instrument 
has  only  been  used,  up  to  November  8,  1991,  to  measure  much  less  intense 
nightglow,  this  was  never  noticed  to  be  a  problem,  and  nothing  yet  is 
conclusive  to  prove  scattering  is  a  problem.  Even  if  scattering  is  the  culprit, 
there  is  nothing  practical  that  can  be  done  to  the  data  to  correct  for  it. 
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Chi-square 


CHI-SQUARE  VS  OI 8446 


18.5  19.0  19.5  20.0  20.5  21.0  21.5 

Local  Time 

Figure  5.1  Plot  of  c2  and  OI  8446  versus  local  time  (EST). 


The  N2  IRA  model  is  another  possible  source  of  improvement.  As 


seen  in  Figure  4.5,  the  modelled  spectrum  does  not  fit  the  experimental 
spectrum  perfectly.  Again  there  was  insufficient  time  to  fully  investigate 
the  model's  output,  so  it  was  necessary  to  use  it  as  is. 


Intensity  (R) 


APPENDIX 


SOURCE  CODE 

The  source  code  for  the  analysis  is  listed  here  for  reference.  The  code 
is  definitely  research-grade  and  was  not  written  specifically  for  optimum 
performance.  Instead  it  was  written  so  that  molecular  and  atomic 
emission  parameters  could  be  added  and  changed  without  a  great  deal  of 
code  manipulation.  Before  listing  the  code,  a  few  helpful  hints  will  be 
passed  along  so  that  someone  else  may  benefit. 

Convergence.  The  first  problem  I  ran  into  was  that  1  could  not  get 
the  chi-square  fitting  routine,  SPEC_FIT,  to  converge:  the  value  of  chi- 
squared  just  got  worse  and  the  recovered  parameters  were  useless.  The 
problem  was  that  when  trying  to  analyze  for  both  band  brightness  am 
temperature  of  a  particular  molecular  species  using  band  models,  you  need 
to  be  in  the  vicinity  of  the  minimum  of  parameter  space.  Both  temperature 
and  brightness  cannot  be  arbitrarily  set.  I  got  around  this  problem  by 
"hard-wiring"  in  reasonable  temperatures  so  that  accurate  brightnesses 
could  be  attained.  Upon  saving  these  parameters  to  a  file,  I  then  adjusted 
the  program  to  read  in  the  file  of  parameters,  and  also  allowed  for 
temperatures  to  then  be  analyzed  along  with  the  other  parameters 
Convergence  was  then  not  a  problem. 
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Data  Spikes.  Another  problem  was  that  the  raw  N1RS  data  had 
"spikes"  in  it.  These  spikes  were  apparently  from  very  energetic  photons 
which  were  able  to  by-pass  the  baffles  and  made  their  way  to  the  CCD. 
Though  only  about  10%  of  the  frames  had  spikes  in  them,  it  did  affect  some 
of  the  fitting  quite  a  bit,  and  the  spikes  seemed  to  be  more  common  during 
the  brighter  periods.  The  spikes  affected  not  only  the  channel  where  the 
photon  was  absorbed,  but  also  about  2  channels  to  either  side  of  the  primary 
channel,  and  many  times  there  were  multiple  spikes  in  a  single  frame  of 
data. 

To  minimize  the  affects  of  the  spikes  on  the  fitting,  any  number  of 
routines  can  be  applied.  I  chose  to  use  a  simple  averaging  scheme.  Once 
the  data  file  was  read  and  returned  to  the  main  program,  I  inserted  a 
series  of  loops  to  check  each  frame  and  channel  for  spikes:  anything 
greater  than  100  counts  over  the  channel's  reading  at  3  frames  to  either 
side  of  the  channel  being  checked,  then  the  5  data  points  within  that  range 
would  be  replaced  by  a  time  average  of  each  channel.  This  time  average 
was  obtained  by  adding  the  preceding  and  succeeding  frame's  reading  for 
that  particular  channel,  and  dividing  by  2.  This  smoothing  process  seemed 
to  work  quite  well.  For  the  worst  fitting  frame  before  the  process  was  added 
had  a  chi-squared  value  of  about  75.  After  the  routine  was  added,  chi- 
squared  was  brought  down  to  under  30,  which  was  close  to  values  obtained 
from  nearby  frames  which  had  no  spikes  in  them. 

Negative  Parameters.  A  third  problem,  though  a  very  minor  one, 
was  that  the  least  squares  fitting  routine  had  no  concept  of  reason-  often 
parameters  (temperatures  and  brightnesses)  would  optimize  to  negative 
values.  I  suppose  negative  brightness  may  imply  absorption,  but  negative 
temperatures  are  difficult  to  grasp.  Though  this  generally  occurred  for 
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only  the  weaker  emissions  and  did  not  affect  the  fit  greatly,  it  did  allow  the 
program  to  "cheat".  This  was  gotten  around  by  resetting  negative 
parameters  to  an  arbitrary,  low  number.  The  value  of  1.0  was  used  in  this 
study.  This  prevented  the  program  from  "cheating"  by  using  better  fitting 
negative  parameters  and  did  not  cause  any  known  problems  on  succeeding 
iterations  when  done  properly. 


c 

C  PROGRAM:  SPEC  ANAL. FOR 
C 
C 
C 
C 
C 
C 
C 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 


PURPOSE:  Main  driver  for  die  SPEC_FIT  subroutine.  First  calls  GETJDATA  subroutine  to 
retrieve  die  NIRS  measurement  data.  It  then  sends  the  data  for  spectral  lilting  to 
SPEC_FIT.  The  fitting  subroutine  uses  a  non-linear  least  squares  iterative  technique  to 
find  the  best  fit  of  die  fitting  parameters  specified  by  the  user  to  the  data.  Recovered 
parameters  are  atomic  emission  line  intensities,  molecular  band  temperature,  :uid  band 
brightness.  Spectral  riuige:  84(X)-X900  A. 

EXTERNAL  SUBROUTINES  CALLED:  MAT1NV,  ATMOS.CROSS,  ( HLSPEC. 
OH_SPECl,  N21P 

MAIN  VARIABLES  USED  IN  PROGRAM: 

NCH:  number  of  channels  of  die  NIRS  instrument 
NPRM:  number  of  panuneters  to  be  recovered  from  die  data 
N:  number  of  atomic  lines  analyzed 
D'f:  integration  period 

MXFRM:  maximum  number  of  friunes  to  be  read  (for  dimensioning  only) 

WAVE:  reference  wavelength.  Determines  where  the  spectrum  hegms  (channel  1 
high  end  since  wavelengUi  decreases  as  the  chtuuiel  number  increases 
CNT:  array  containing  die  number  of  counts  (corrected)  amuiged  by  diaiviel,  frtune  no 
EMI:  synthetic  spectrum  counts  which  best  fas  the  data  based  on  specified  parameters: 

output  from  SPEC_FJT  to  compare  with  data  (COUNT). 

COUNT:  ;ui  array  of  the  current  frames'  data  from  CNT  to  be  sent  to  SPHCJTI . 

A:  array  of  recovered  panuneters.  Dimensioned  to  NPRM. 

SIGMA:  array  of  errors  of  die  A  pttrtuneiers.  Dimensioned  to  NPRM. 

LINES:  turay  of  atomic  lines  to  be  analyzed. 


at  the 


Designed  by: 
Written  by: 


Thomas  .1.  Smith 
Thomas  J.  Smith 


15  January  1992 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCOCCCCC CGCCCCCCCCCCCC 


PARAMETER  (NCH=384.  NPRM=3l,  N=ll,  DT=3().().  MXFRM=75(» 
REAL* 8  WAVE 

REAL  CN'r(NCII+l.MXI'RM).EMI(NCU).TIME(MXFRM).CCD(MXIRM) 
REAL  COUNT(NCH),DK(MXFRM),A(NPRM).SIGMA(NPRM).LINl-:S 
INTEGER  SPELG(MXFRM) 

COMMON  /OH/ID_OII.DARK,LINE.S(N) 

COMMON  /TRANS/WAVE 
DATA  WAVE/  X906.249D0/ 

DATA  LINES/  8446.50.  8567.74,  8594  (H).  8629.24,  8655.89. 

>  8680.28.  8683.40,  8703.25.  871 1.70.  8718.83,  8828.89/ 


C  Allow  for  different  ( )11  molecular  data  to  be  used. 


C  WRrrE(6,1020) 

C  1020  FORM AT(  1  X.’ENTER  THE  SOURCE  OE  OH  TRANSITION  PROBABILITY:./. 
C  &  1X,'0:  MIES;  ETURNBULL  AND  LOWE:  2:  LANGHOFF:  \S) 


n  n 


C  ACCEPT  *,  lD_OH 
ID_OH  =  0 


C  Get  the  data  to  analyze 

CALL  GET_DATA(DT,  CNT.Nl  IMFRAM.T1ME.DK.CCD) 

PRINT*,'  FRAME  TIME  CHI  SQ  CCD  TEMP  U  LOOPS  ' 

C  Now  check  data  lor  spikes  by  finding  an  avg  count  from  3  channels  either  side,  and  avg  of  l 

C  channel  either  side  of  each  point  If  the  data  is  greater  than  100  counts  (RUM)  more  than  the  2 

C  averages  then  the  data  point  is  assigned  the  value  of  die  time  average  (AVI'RAM). 


RLIM=100. 

DOI=2,NUMFRAM-l 
DO  J=3,NCH-3 

AVCHAN  =  (CNT(J+3,I)+CNT(J-3.I))/2. 

AVFRAM  =  (CNT(J,I+l)+CNT(J.I-l))/2. 

IF  ((CNT(.I.I).(iT.AVCHAN+RLIM).AND. 

&  (CNT(J.I).GT.AVFRAM+RLIM))  THEN 

IF  ((CNT(J.I).GT.CNT(J+LI)+RLIM).AND. 

&  (CNT(J,l).GT.CNT(J-l,I)+RLIM»  THEN 

CNT(J-2,I)  =  ( CNT(  J  -2.1+ 1  )+CNT(J-2.1  - 1  ))/2. 

CNT(J-I.I)  =  (CNT(J- 1 .1+1  )+CNT(J- 1 ,1- 1  ))/2. 

CNT(J.I)  =  AVI'RAM 

CNT(J+1,I)  =  (CNT(J+ 1 ,1+ 1  )+CN  1  (J+ 1 .!- 1  ))/2. 

CNT(J+2,I)  =  (CNT(J+2,I-t  1  )+CNT(J+2,I- 1  ))/2. 

SPFLG(I)  =  1 
ENDIF 
ENDIF 
ENDDO 
ENDDO 

Loop  over  die  intinber  of  frames  of  data  (Nl  IMFR  AM)  for  antdysis.  Reassign  the  array  ol 
intensities  CN  F  to  the  vector  COUNT  then  pass  it  to  SPEC_FIT. 

5  DO  1(X)  l=LNUMFRAM 
DO  10  J=1,NCH 
10  COUNTED  =  CNT(J,I) 

DARK  =  DK(I) 

CALL  S  PECJ  TT(  COUNT.  DT,  A.SIGMA.EMI.NLOOPS.XSI.IERR) 

IF  (I.EQ.I)  THEN 

OPEN(1  FN1T=2,FIL  E='[T. SMITH.  AURORA.  DA' TA]3 1  PRMSPI-C.DAT. 

&  ST  ATI  IS='NEW') 

OPEN(UNJT=3,FILE=  ! TSMITH  AI  IRORA.DATAJ31PRM.DAT, 

&  ST  ATI  !S=‘NEW') 

ENDIF 

IF  (IERR.EQ.O)  THEN 

WRITE(6,’(10X,I5,3F10.3,2I8),)I,TIME(I),XSI,CCD(I).NL(X)PS,SPn,G(I) 

WRJTE(3,,(F8.3.F7.2,F6.L2I3,10F8.I)')TIME(l),XSI.OCD(I), 

&  SPFLG(I),Nli)OPS.(A(J).J=I,10) 

WRITE(3,'(27X,1 1 1-8. 1 )')( A(J>,J=  1  l.NPRM) 

ENDIF 

100  CONTINUE 
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noorinnnn 


CLOSE  (2) 
CLOSE  (3) 
STOP 
END 


ccccccccccccceccccc.  SUBROUTINE  GETJMTA  cccccccccccccccc.ccccc.cc 

SUBROUTINE  GET_DATA(DT,  CNT.Nl IMERAM.TIM1%DARK.CCD) 


PURPOSE:  Reads  in  die  NIRS  speciromeier  data  file  by  tune  (triune)  then  corrects  die  data  tor 
dark  counts  and  instrument  bias  (mad  error). 


MAIN  VARIABLES: 

DARK:  array  of  dark  counts  measured  by  die  (covered)  TKSdi  channel  during  each  frame 
CCD:  array  of  CCD  temperatures  calculated  from  die  dark  counts 

CNT:  array  of  counts  measured,  corrected  for  dark  counts  and  bias.  Arranged  by  channel,  triune 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


PARAMETER  (NCI  1=3X4.  MXERM=75(» 

REAL  CNT(NCH+ 1  ,MXERM),CO(NCI  1+ 1 ),('  1  (NCI  1+ 1  ).I  KNCi  1+ 1 ) 

REAL  TIME(MXERM),DARK(MXERM).C(  H  JNT(NCl 1+ 1  .MXERM).CCD(MX'-RM) 
INTEGER*4  STATUS,LIBSGET_U  IN 
INTEGER  DARKJJNIT.BIAS _l I) JIT. IK )l IR.MIN 


C  First  open  instrument  calibration  and  measurement  Jala  tiles 


STATUS  =  LIBSGET  LUN(DARK_UNIT) 

OPEN  (UNIT  =  DARKJJNIT, 

&  N  AME='(TSM  ITI I.  AURORA  DATA  [MAI  7. S'  DARK_('OE'E  3H4.DAT. 

&  READONLY,  TYPE='OLD) 

STATUS  =  LIB$GET_LUN(BIAS_UNIT) 

OPEN  (UNIT  =  BIAS_I  JN1T. 

&  NAME='[TSMITH.AURORA.DATA]BIAS  384.DAT, 

&  READONLY.  TYPE='OLD') 

OPEN  (UNIT  =  1, 

&  NAME=’(TSM  [TIE  AURORA. DATAJ.SPECT.DAT, 

&  NAME=’SPRI.  ,I::HRDSSCR:(  YEEJNOV08yi.DAT. 

&  READONLY,  TYPE=  ()LD') 

PRINT  ' 

PRINT  READING  INST  RUMENT  DATA...' 


C  Read  the  dark  counts  and  bias. 


DO  K  =  l.NCH+1 

READ(DARK_UN IT,  *)  C  J IM I P.C  I  ( K ). H( K ) 
READ(BIASJ  JNIT.*)  KK.('0(K).C()_ERR 
ENDDO 


NUMFRAM  =  « 

TIMEJJLD  =  OW.O 
INDEX  =  0 

PRINT  V  READING  DATA  ETI.E 
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C  Loop  to  read  in  ine:uisuremen(  data,  check  tor  had  tnune  (('( )l  JNT=163X3) 

DO  900  I  =  1.MXFRM 
NUMFRAM  =  NUMFRAM  +  1 
READ(1.*,END=950)  lA.IB.IC.ITlME.ICH.CCT.N.TP. 

&  (COUNT(j,NUMFRAM)J=LNCH+I) 

IF  (ITIME.LT.  12(H))  ITIME  =  ITIME  +  24(H) 

IF  (ITIME.LT.  1X00)  THEN 
NUMFRAM  =  NUMFRAM  -  1 
GO  TO  900 
ENDIF 


DO  J=  l.NCH+l 

IF  (COUNT(J,NUMFRAM).EQ.  16383)  THEN 
NUMFRAM  =  NUMFRAM  -  1 
GO  TO  900 
FNDIF 
ENDDO 

HOUR  =  ITIME/ 100 
MIN  =  ITIME  -  HOUR*  100 
MIN  =  HOUR*60  +  MIN 
TIME(NIJMFRAM)  =  MIN/60.0 
IF  (TIME(NUMFRAM).EQ.TIME_OLD)  THEN 
TIME(NUMFRAM)  =  TIME(NUMFRAM)  +  INDEX  *0.008333 
INDEX  =  INDEX  +  1 
ELSE 
INDEX  =  1 

TIME  OLD  =  TIME(NUMFRAM) 

ENDIF 

900  CONTINUE 

950  NUMFRAM  =  NUMFRaM  -  1 
CLOSE  (1) 

CLOSE  (DARKJINIT) 

CLOSE  (BIAS_UNIT) 

PRINT  *.  ’  TOTAL  FRAMES  READ  =\Nl  I  MLR  AM 


C  Subtract  out  the  instrument  bias,  then  calculate  CCD  temp  and  DARK  count,  then  subtract  out  the 
C  dark  count. 


DO  I  =  1, NUMFRAM 
DO  J=  l.NCH+l 
CNT(J.I)  =  COUNT(J.I)-  OK385) 
ENDDO 


CCD(I)  =  -  50.0  -  H(3X.5)*LOG(CNT(3X5,I)/(Cl(3X5)*DT)) 
DARK(I)  =  DT*C  1  ( 1  )*EXP(-(CCD(l)+50.0)/I  K 1 )) 


DO. I  =  l.NCH 

CNTU.I)  =  CNT(J.l)-DT*Cl(J)*EXP(-(CCD(I)+5().())/H(j)) 
ENDDO 
ENDDO 
RETURN 
END 
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S  UB  ROUTINE  SPEC  JIT  CXXXXX X’(X XX XX X ( X XX XX X ' 


si  iBROirriNU  SPECJH(Y.1>T,  A.SIGMA.E.NI.GOPS.XSI.IERR) 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

n 

V- 

c 

c 

c 

c 

cccr 


PURPOSE: 

Non-linc;ir  least  squares  lilting  of  data  to  derive  par.uneters  such  brightness  and  hand  temperature. 
Uses  an  iterative  technique,  halting  iteration  when  convergence  ol  each  parameter  ol  A  is  within 
0.5%  of  previous  value.  Obtains  the  synthetic  spectrum  ol  partuneters  by  sending  them  to 
EMLSPEC  which  uses  NlRS_CONVOEVER  to  convert  synthetic  spectrum  to  what  N1RS  would 
measure  at  each  channel. 

MAIN  VARIABLES: 

Y:  array  of  measurements  of  dimension  NCI  I 

NPRM:  number  of  pturuneters  to  be  derived-  determines  matrix  dun’s 

NCH:  number  of  channels  in  spectrum-  determines  summation  size 
A:  array  of  derived  partuneters  to  be  output  when  iteration  is  complete 
SIGMA:  array  of  situulard  deviations  of  A  partuneters 

F:  array  of  emissions  predicted  by  EMUS  PEC  based  on  inputted  parameters,  i.e.  ft  A) 

DFDX:  array  of  eenlered-differeuee  derivatives 

NOTE:  tile  resolution  used  in  the  subroutine  is  set  by  the  parameter  RE-S.  ( )nce  convergence 
occurs  to  less  than  RES  for  till  partuneters,  the  loop  halts  and  A  is  returned. 
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCi  cccccccccccccn  c 


PARAMETER  (RES=1.<).  NPRM=31.  N=l  1.  NCIU0X4.  MAXL(X)PS=10) 


REALMS  DFDX(NPRM,NCII).B(NPRM),ARRAY(NPRM,NPRM) 
REAL  Y(NCH),DEL(NPkM),A(NPRM).SIGMA(NPRM),E(NCH) 
REAL  CNT2(NCHLX(NPRM).YI(NC11).Y2(NCH).1)X(NPRM).LINES 
INTEGER  INDX(NPRM) 

CHARACTER*!  CONV 

COMMON  /Ol f/ID  Ol  I.l)ARK.LINliS(N) 

LOGICAL*  1  ElRST/.Tkl  IE./ 

1ERR=() 

NCOL  =  1 
NLOOPS  =  0 


C  Initialize  partuneters  for  a  reasonable  first  guess 
IF  (FIRST)  THEN 

OPEN (l JNn'=4,FIEE>'[TSMITll. AURORA  DATA|31PRM  DAT'.STATl  LS=  C  >1.1)  ) 
READ(4,*)TIME.CHI,CCI)TEMP,L( )( IP 
READ(4,*)(K,A(I),ERR.1=1,31) 

CLOSEC4) 

FIRST  =  .EAI.SE. 

ENDIF 

CALL  EM1_SPEC(I)T,N(  OI..A.O,  id 

C  Estimate  die  accuracy  of  the  observed  spectra  (photon  statistics  plus  the  bias  statistics  ( 1 .2)). 
C  CNT2  is  the  vtuiance  of  each  channel. 

1202  1=  l.NCII 

2  CNT2(I)  =  (Y(I)+I)ARK)/15.6+  1.2 
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CCCCCCCCCt 'CCCCCCGC  MAIN  LOOP  CCCCCCGCCCOCCCt '( ( X '( VC 
5  NLOOPS  =  NI.OOPS  +  I 

C  Calculate  the  derivative  array  DX  usiti”  centered  dilterence  method. 


LX)  140  l=l.NPRM 
DO  10  J  =  1,NPRM 
10  DX(J)=0,0D0 

DELTA = RES 
DO  100  J=  1,2 
DELTA  =  -DELTA 
DX(I)  =  A(t)  +  DELTA 


0  Set  up  derivative  arrays  so  lltat  each  element  except  lor  needed  elements  is  zero. 

C  Dx  is  +/-  1  R  or  I  K. 


& 

& 


16 

17 

20 

21 

24 

25 


28 

20 


32 


C 

O 

33 

80 

100 


GO  T( )( 1 2, 14, 80, 80,80.80,80. 80,80.80, 
80,80.80,16,16. 16,17,20,28.28, 

2 1,24,24.25,20,32, 33,33, 33, 16,80),  1 
DX(2)  =  A(2) 


GO  TO  80 
DX(  i )  =  A(  1) 

GO  TO  80 
DX(N+6)  =  A(N+6) 
GO  TO  80 
DX(N+3)  =  A(N+3) 
DX(N+4)  =  A(N+4) 
DX(N+5)  =  A(N+5) 
DX(N+2<»=  A(N+20) 
GOTO  80 

DX(N+10)  =  A(N+10) 
GO  TO  80 
DX(N+7)  =  A(N+7) 
GO  TO  80 

DXCN+I3)  =  A(N+ 13) 
GOTO  80 

DX(N+1 1 )  =  A(N+1 1) 
DX(N+12)  =  A(N+ 12) 
GO  TO  80 

DX(N+14)  =  A(N-t-14) 
GO  TO  80 
l)X(N+8)  =  A(N+8) 
DX(N+0)  =  AlN+0) 
GO  TO  80 

DX(N+16)  =  A(N+16) 
I)X(N+17)  =  A(N+17) 
DX(N+ 18)  =  A(N+18) 
DX(N+1<))=  A(N+ld) 
DX(N+20)  =  A(N+20) 
GO  TO  80 

DX(N+15)  =  A(N+ 15) 


!  CONTINUUM  1 

!  CONTINUUM  2 

!  N21P  BANDS 

!  (1.0) 

!  (2.1) 

!  (3.2) 

!  (4,3) 

!  02  A TM(  )S  BANDS 

!  (0.1) 

!  Oil  BANDS 

!  (6.2) 

!  (7,3) 

!  02  (1,2)&(2,3)  TEMP 

’  (1.2) 

!  (2,3) 

!  N21RA  BANDS 


IF(J.EQ.l)  GALL  EMI_SPEC(DT,NCOf „DX.I.  Yl) 
IE(.I.EQ,2)  CALL  EMLSPEC(DT,NC<>L.DX,1,  Y2) 
CONTINUE 


DO  130  1=1, NCI  I 

130  D1'DX(U)=(Y2{.I)-Y1(J))/(2.0D0*RES) 
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140  CONTINl  Hi 


DO  170  1=1.  NORM 
DO  160  J=I,NPRM 
ARRAY(U)=O.ODO 
IX)  150  K=l,NOH 

150  ARRAY(U)=ARRAY(U)+ 

&  DFDX(LK)*DEDX(J.K)/CNT2(K) 

160  CONTINUE 

B(1)=0.0D0 
DO  165  L=1.NCII 

165  B(I)=B(1)+(Y(L)-E(U)+DFDX(1.L)/CNT2(1.) 

170  CONTINUE 

C  Invert  ARRAY  matrix,  tlien  solve  tor  X=ARRAY*B 

CALL  MATINV( ARRAY, NORM, Nl’RM) 

DO  550  1=1, NORM 
XU)  =  0.01)0 
DO  500  J=l, NORM 

500  X(l)  =  X(l)  +  ARRAY(U)*B(.T) 

DEL(I)  =  X(I) 

A  (I)  =  X(I)  +  Ad) 

550  CONTINUE 

C  Calculate  error  of  each  parameter 

DO  580  I=1,NPRM 
IF(ARRAYU.I)  OR.0.01 )()) 

&  SICiMA(l)=SQRT(ARRAY(l.D) 

IF(A(1).LT.  1  .OE-5)  THEN 
DEL(I)  =  0.0 
A(I)=  1 .0 
END1F 

580  CONTINUE 

C  Calculate  error- weighted  Chi-squared 

CALL  EMI_SPEC(DT,NCOL.A,0, 1) 

XSI  =  0.0 
DO  I  =  1,  NCH 

XSI  =  XSI  +  (Y(I)-RO)**2/CNT2(I) 

ENDDO 

XSI  =  XSI/(N011-NPRM) 

C  Now  test  each  parameter  tor  convergence  using  0.05 %  as  criterion.  Repeat  loop  il  convergence 

C  has  not  occurred  (lor  parameters  w /  magnitude  >  1  (X))  unless  MAXLOOPS  iterations  have 

C  occurred. 

CONV  =  *Y’ 

DO  600  I=1,NPRM 

600  IF  (ABS(DEL(1)/A(D)  ( IE. 0.005  AND. 

&  A(I).CE.  100.)  C( )NV='N‘ 

IF  (NLOOPS.LT.MAXLOOPS  .AND. 

&  CONV.EQ.’N')  CiOK)  5 

999  RETURN 
END 
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'  For  use  in  convergence  test 
!  SOLVE  lor  u(i)  for  use  in  next  iteration 


SUBROUTINE  EMI  SPEC  ({'('{<<{  CCCCCt X  C( XXXX  C 


SUBROUTINE  !;M1_SPEC(DT.NC(  )E,A.1ST.  EMI) 

C  PURPOSE: 

C  To  calculate  an  emission  spectrum  based  on  A  parameters  that  c;ui  be  used  to  compare  with  a 

C  measurement  from  the  NIKS  instrument.  Emission  spectrum  is  sent  to  NIRS_CONVOEVER 

C  subroutine  to  convert  Ihousttnds  of  atomic  and  molecuittr  lines  to  a  channel-hv-ch;utnei  signal, 

C 

C  MAIN  VARIABLES: 

C  NCOL:  number  of  pixels  per  column,  or  channel,  set  on  the  (XT)  at  measurement  tune. 

C  W:  array  of  wavelengths  returned  from  b;uid  models 

C  WET:  array  o!  wavelenghts  of  all  emission  lines 

C  S:  array  of  strengths  relumed  from  ha; id  models 
C  STR:  turay  of  strengths  of  all  emission  lines 

C  Nl  IM 1 :  number  of  lines  returned  from  a  haul  model 

C  Nl  !M:  total  number  of  lines  to  be  convolved 


PARAMETER  (NC!!=3X4.  NPRM=31.  N=1 1.  MAXEN=1(XXH)) 
REAL  LINES,  A(NPRM),EMI(NC11).WK2<K)).W2(2<X)) 

REAL*  8  WAVE.W(  15(M)),S(  I5(X)),WLT(M  AXEN).STR(MAXEN) 
COMM( )N  /< )1 1/11  ),1)ARK,1 .INES(N) 


C  Call  the  appropriate  band  model(s)  given  by  Uie  1ST  indicator 


ISTEP=IST+I 

NUMTOT=() 

GO  TOO  1,‘)().‘)(),()().‘)(),9().‘)(),()().V().<)(). 

&  ‘)().<)().<K),y(),  1 1,12.13,11.21.22. 

&  23, 2 1,3 1,32, 3 1.22, 4 1,4 1,42.43. 

&  14. 52), INTER 

1 1  CALL  N2 1  P(  1  ,(),A(N+6),W„S,N(  IM  1 ) 

CALL  ADDLIN(N+3, 1.N1  !M  1  .Nl  JMTOT.W.S.A.  WLT.STR) 
IF  (ISTEP.NE.  1  .AND.  ISTEP.NE.N+7)  GO  TO  90 

12  CALC  N2 1  P(2, 1  ,A(N+6),W.S,NUM  1 ) 

CALL  ADDLIN(N+4. l.Nl  IM  1  .NUMTOT.W.S.A.  WLT.STR) 
IF  (ISTEP.NE.  1  AND.  ISTEP.NE.N+7)  GO  TO  90 

13  CALL  N21P(3,2.A(N+6).W,S.N1  IM  1 ) 

CALL  ADDLIN(N+5. 1  .Nl IM  1  .Nl  IMT( )T.W.S.A,  WLT.STR) 
IF  (ISTEP.NE.  1  .AND.  ISTEP.NE  N+7)  GO  TO  90 

14  CALL  N21P(4.3.A(N+6).W.S.NUM1) 

CALL  ADDI.IN(N+  PLENUM  I, NUMTOT  .W.S.A.  WET  , S  I  R) 
IE  (ISTEP.NE.  I)  GO  TO  <>(> 

2 1  CALL  AITVK )S_GR( )SS((),1 , 1 , A(N+  10),W.S,W  1  .W2.NUM1 ) 
CALL  ADDLIN(N+7, l.NUM  1  .Nl  JMTOT.W.S.A.  WLT.STR) 
IF  (ISTEP.NE.!)  GO  TO  DO 

22  CALL  ATM( )S_CR( )SS(  1 .2, 1 , A(N+  14).W,S.W  1  .W2.N1  IM  1 ) 
CALL  ADDL1N(N+S,  l.NUM  l.Nl  IMTOT  .W.S. A.  WLT.ST  R) 

IF  (ISTEP.NE.  1  .AND.  ISTEP.NE.N+1 5)  GO  TO  90 

23  CALL  ATMOS  CROSS(2.3.I.A(N+I4).W.S,WI,W?,NOMI) 
CALL  ADDLIN(N+9.  l.NUM  1. NUMTOT.W.S.A.  WLT.STR) 

IE  (ISTEP.NE.  I)  GO  TO  VO 

31  CALL  Oli_.SPE(’(ft.2.A(N+l  1).A(N+E3).I1).  W.S.NUM1) 

CALI.  ADDEIN(N+1 1.2.NUM1. NUMTOT.W.S.A.  WLT.STR) 
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IF  (IS'lEP.NLl  AND.  1STEP.NE.N+14)  GOTO  oo 
32  CALL 01L.SPHCl(7,3.A(Ni  12).AtN+I3).ID.  W.S.NUMI) 

CALL  ADDL1N(N+  12.2.NI  IM 1  .Nl  !MT(  )T.W.S, A,  WLT.STR) 
IF  (ISTEP.NE!)  GO  TO  00 


41  CALL  BANDS  YN(6,A(N+lf)),A(N+ 1 5).  1 1 .6.W.S.NI IM  1 ) 

CALL  ADDI ,IN(N+16,3,NI  IMl.NlFMTf VL.W.S.A,  WLT.STR) 

IF  (ISTHP.NE.  1  .AND.  ISTEP.NE.N+lfi)  CiO  TO  40 

42  CALL  BANDSYN(6,A0N+17).A(N+15),6,2.W,S.NUM1> 

CALL  ADDLIN(N+ 17,3, NUM  l, NUMTO  L.W.S.A,  WLT.STP.) 

IF  (ISTEP.NH.l  .AND.  ISTEP.NE. N+ 16)  GO  TO  00 

43  CALL  BANDS  YN(6,A(N+ IX).  A(N+ 15),  12,7,  W.S.NUM  1 ) 

CALL  ADDLIN(N+ 1 8,3, NI  IM1.NIJM TOT.W.S.A,  WLT.S  1  R) 

IF  (ISTEP.NE.  1)  C>()  TO  00 

51  CALL  BANDSYN(5,A(N+20),500  (),4,2, W.S.NUMI) 

CALL  ADDLIN(N+20,3,Nl  IM  1  .Nl  !MT< )  l  .W.S.A.  WI .1  .S  I  R ) 

IF  (ISTEP.NE.  I )  (i()  T( )  00 

52  CALL  BANDSYN(5,A(N+2 1.500.0. 5. 3.W.S.NUM1) 

CALL  ADDLIN(N+20,3.N(  IM  I  .NI  1MTX )  L.W.S.A,  WI  T  ST  R) 

00  CONTINUE 

C  Now  add  (he  atomic  linos  (o  die  WLT  and  S  LR  arrays  diet)  convolve  die  arrays  into  a  channcl-hy- 
C  channel  signal. 


DO  1=1, N 

WL.TCNl  JMTOT+1)  =  LINES!  1) 

STR(NUMT()T+1)  =  A(I+2) 

ENDDO 

NUM  =  NUMTOT+  N 

CALL  NIRS_(’()NV( )LVER(WI .T.STR.Nl  IM.Al  1  ).A(2).NC(  )L.DT.  EMI) 

RETURN 

END 


CCCCCCCCCCCCCCCf  XX XX  X  ’  SI  ro R< )l !  LINE  AI )DI  .IN  C  X  X X X  X XXX  X X X X X X X X  XX X XX’ 
SUBROUTINE  ADDLINdIX )l IFLAC  i.NI IM  i  .Nl  JMT( )  L.W.S.A.  WLT.STR ) 

C  PURPOSE: 

C  Adds  emission  lines  from  the  W  and  S  arrays  to  die  WI.T  and  S  LR  arrays  respectively,  but  first 

C  convert  die  different  strengths  given  by  the  dillerent  band  models  to  a  band  brightness 

C  normalization  il  necessary. 

C 

CCCCCCCCCCCCCCCCCCCCCCIX  (X XXX  (XX  CCCIXXX  (  (XX  CCCOCI X  (XXXXXX X  cc 


PARAMETER  (NPRM=31.  MAXI ,N=  10000) 

INTEGER  OI  ILL  AG 
REAL  A(NPRM) 

REAL*X  W(  1 500). S(  1 500), WI  .KMAXLNLS  LR(MAXLN) 

GOTO  (10,2(),3(».OI II  LAC. 

10  DO  15  1=1, NUM  1 

WLT(NUMT()T+I)  =  W(l) 

15  STR(NUMT(  )T+I)  =  S(I)*A(II) 

CiO  TO  40 

20  DO  25  1=1, NUM  1 

WLT(Nl IMTOT+I)  =  1.01  >3/(1  000275D-5*W(1)) 
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25 


ST1<<NUMT<  >1+1  >  =  S(D 
GO  TO  40 
30  DO  35  1=1.NUMJ 

Wi;r<NlIMTOT+I)  =  \V(l) 

35  STR(NUMTOT+l)  =  S(l> 
40  CONTINUE 

NUMTOT  =  NIIMTOT  +  NUM1 

RETURN 

END 


subroutine  nirs.gonvoever  oca  xx  x  < x  t <  <  <  <  <  < xxxxx 


SUBROUTINE  NikSJONVOI  VERtWAVEl.D.NUM CONT.GON2.NGO!  .01.  SIG) 

(’  PURPOSE: 

C'  lakes  input  ol  wavelengths,  dictr  associated  line  strengths,  and  high  and  low  continuum  .and 

('  convolves  the  lines  into  a  chanucThy-clinnuel  signal  that  would  lx-  measured  h>  the  NIPS 

(’  instalment 

G 

('  DESIGNED  BY  JENG-I1WA  Yi  l 
G  WRITTEN  BY  1ENG-IIWA  YEP. 

C 

C  1)ATE=  ( XT  2.  POO 

CCCCCGCCC X '( X 'GC '( '(  X  X  X  X  X  X  X  X  X  X  X  XXX  '(  X  X  XXX  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  (  (  <  <  <  <  X  X  ( 

PARAMETER  (MAX1.N=  10000.  NCH=3X4> 

REAE*X  WAVT..VVAVEE<MAXEXj.B(MA\T  .NT 
DIMENSION  SPEC(5i>0>.SPEC  P(500) 

DIMENSION  S IGN Al .( NCI  I  ).S  1G( NCI  I ) 

DIMENSION  TR(27).SEN(N('1I) 

LOGICAL*!  EIR.STJ’AEI./  TKI  IE./ 

DATA  ADJINlT/I5/>/ 

C(.)MM<  )N  /trans/wave; 


ie(First_gau.)'hu:n 

OPENaiNrr=II.l;IU>  ['l. SMITH. AURORA  DATAjMAHS  INST  01  <>|  DAE. 
&  READONLY  .SHARED.STATl  IS=X  )l  TV) 

OPEN(l  INIT=12,iTI.i-=’(  ISM.'I  II.AI  IRORA  DATA  [MATES  SPINS  01  01. DAT. 
&  RHADONEY.SIIARII  ).S  1  AT  l  >S=X)I  IV) 

DO  .1  =  1,27 
READOI.O  Hid ) 

ENDDO 

DOE  =  INCH 

READ02,*)  WAN'I  l  I !  SI  Ntl  .i.ST.NS  ERR 
ENDDO 

DW_DP  =  14035 

dw  =  dw_dp*eoe:x/<wave;-ioo.o)**2 

TOT  =  0.0 
DO  I  =  1,27 

TOT  =  TOT  +  TR(I)*DW  DP 
ENDDO 


PIRST_CAPP  =  1A!  SI 
PND1P 


S!  A )I*i:  =  (CON2-C '( )N  1  )/i  I .< ) A  I  (NTH) 
DC)  I.  =  1.500 

SPPX'(P)  =  {('( )N1+  (I. •'>((»'.' ’  •  ,  r)-D\V 
KNDDO 


DO  I.  =  l.NDM 

II  =  (WAVT  -  '"AVi:i.(i.)>/D\\  DP  *  50 
IP  (II  LPT  AND.  1 1  OT.0)  SPPOIH  -  SPJ-(’III)  +  Dfl  ) 
PNDDC ) 

1X>  I.  -  l.NC'II 
SICiNAI.(L)  =  0.0 
SP1X'_P(L)  =  0.0 
PNDDO 

DO  I  =  l.NC'II 
II  =  I  +  50 
DO. I  =  1,27 

SPP('_P(I)  =  SPIT  Pd)  +  TKi  I 1 * N PI •  ( '( 11-14  + 1 1 
PNDDO 
PNDDO 


DO  I.  =  l.NC’II 
II  =  <P- 1  )*N( '( >1.  +  i 
DO  I  =  II.  II+NC'OP-I 
SKiNAP(P)  =  SIC  iNAI.d  )  +  S  PI  A  hit 
PNDDO 
PNDDO 
DO  N  =  l.NC'II 

S 1 C  i(  N  )= S  PN I N )  *  S I O  N  A !.( N )  •  D  P(  1 0 1  •  N  <  OP  I 
PNDDO 
RP'II  !RN 
PND 


!>} 
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